$/£***- 

/Ay- 3^ - e/2 - — 



The Ohio State University 


ELECTROMAGNETIC SCATTERING FROM A CLASS OF OPEN-ENDED 
WAVEGUIDE DISCONTINUITIES 
By 

A. Altintas 
P.H. Pathak 
W.D. Burnside 


The Ohio State University 

ElectroScience Laboratory 

Department of Electrical Engineering 
Columbus, Ohio 43212 


Technical Report 716148-9 
Grant No. NSG 1613 
March 1986 


National Aeronautics and Space Administration 
Langley Research Center 
Hampton, Virginia 23665 


(NAS A -CR- 181271) 
f EC M A CLASS CP 
E1SCC NTINUITI ES 
Avail: N1IS EC 


ELECTBOBAGNEIIC SCAT1EBING 
CEEN-EKDED BAVEGUJDE 

(Chic State triv.) 319 p 
A 14/BF A0 1 CSCL 20N 


H87-27874 


G3/32 


Unclas 

C093199 



NOTICES 


When Government drawings, specifications, or other data are 
used for any purpose other than in connection with a definitely 
related Government procurement operation, the United States 
Government thereby incurs no responsibility nor any obligation 
whatsoever, and the fact that the Government may have formulated, 
furnished, or in any way supplied the said drawings, specifications, 
or other data, is not to be regarded by implication or otherwise as 
in any manner licensing the holder or any other person or corporation, 
or conveying any rights or permission to manufacture, use, or sell 
any patented invention that may in any way be related thereto. 


30272 MQ) t 

REPORT DOCUMENTATION 
PAGE 


i. REPORT HO. 


2 . 


4. Title end Subtitle 

ELECTROMAGNETIC SCATTERING FROM A CLASS OF OPEN-ENDED 
WAVEGUIDE DISCONTINUITIES 


7. AuthorU) 

A. Altintas, P. Pathak, W. Burnside 


Performing Organization Name and Address 

The Ohio State University 
ElectroScience Laboratory 
1320 Kinnear Road 
Columbus, Ohio 43212 


12. Sponsoring Organization Name and Address 


National Aeronautics and Space Administration 
Langley Research Center 
Hampton, Virginia 23665 


3. Recipient's Accession Nt* 


5. Report Date 

March 1986 


B. Performing Organization Rapt. No. 

716148-9 


10. Project /Task /Work Unit No. 


11. Contract(C) or Grant(G) No 
(C) 

(O NSG-1613 

13. Type of Report & Period Covered 

Technical 


M. 


15. Supplementary Notes 


IS. Abstract (Limit: 200 words) 


A relatively simple high frequency analysis of electromagnetic scattering from a 
class of open-ended waveguide discontinuities has been developed. The waveguides are 
composed of perfectly-conducting sections in which the electromagnetic field can be 
written as a sum of waveguide modes. Junctions are formed at the open-end and also 
within interior regions where different sections are joined. The interior modal field 
is expressed in terms of an equivalent set of "modal rays". The reflection and 
transmission properties of each junction are described in terms of a scattering matrix 
which is determined by combining the modal ray picture with high frequency techniques 
such as the geometrical theory of diffraction (GTD) , the equivalent current method 

(ECMV a«d mnHi fi rat i fine of the nhusi ral thar.r'j rtf HiffrarHnn fPTH^ ft now cot nf 

\ /» * “ w . vi.» ‘ r ’ ’ J ” ’ * W ‘ J ~ ’ - . . . . ^ v . . w> , y . ■ ^ / . ' • < 

equivalent currents are employed in this ECM analysis which leads to a simple treatment 
of many types of junction discontinuities. Also, a new procedure is presented to 
improve the efficiency of the aperture integration at the open end which is required in 
the PTD procedure for finding the fields radiated from (or coupled into) the open end. 
Once the scattering matrices are determined via the aforementioned high frequency 
techniques, they are then combined using a self-consistent multiple scattering method to 
obtain the total scattered fields. The accuracy of the present approach has been 
verified by comparison with other available solutions and measurements wherever possible. 


17. Deeument Analysis a. Descriptors 


b. tdantlflers/Open-Ended Tarmt 


c. COS ATI Field/Group 



19. Security Class (This Report) 

Unclassified 


20. Security Class (This FW) 

unclassified 


21. No. of Pages 

301 


22. Price 


(Sea ANSUZ39.1S) 


See Instructions on Reverse 

i 


OPTIONAL FORM 272 (4-77) 
(Formerly NTIS-3S) 
Department of Commerce 




TABLE OF CONTENTS 


PAGE 

LIST OF FIGURES v 

CHAPTER 

I INTRODUCTION 1 

II SELF CONSISTENT MULTIPLE SCATTERING MATRIX FORMULATION 13 

III DEVELOPMENT OF THE ELEMENTS OF THE SCATTERING 

MATRICES IN THE MULTIPLE SCATTERING METHOD 27 

3.1 The Reflection Type Scattering Matrix 27 

3.1.1 Reflection Back Into the Exterior Region 28 

a) Near field scattering by the open end of a 

c pmi -*i nf -j n-i t o circular W3V0QUidGt 34 

b) Near field scattering by the open end 

of a semi -infinite rectangular waveguide: 46 

3.1.2 Modal Reflection from an Interior Discontinuity 50 

a) Modal reflection at the open-end of a 

parallel -pi ate waveguide: 61 

i) TM case: 62 

ii) TE* case: 62 

z 

i ii 


PRECEDING PAGE BLANK NOT FILW=D 



b) Modal reflection from the open-end of a 

rectangular waveguide: 63 

c) Reflection of TE 1(1 mode from an E-plane 

circular bend in a rectangular waveguide: 68 

d) Modal reflection from a junction between two 

linearly tapered waveguides: 73 




case: 

case: 


e) Modal reflection from the junction of two 
sectoral waveguides: 


73 

74 

75 


i) TE X case; 76 

ii) TM X case: 77 

f) Reflection from the open end of a circular 

waveguide: 79 

i) TE case: 79 

ii) TM* case: 31 

g) Modal reflection in an annular waveguide 

terminated by a paral lei -plate waveguide: 85 


The Transmission Type Scattering Matrix 90 

3.2.1. Transmission Between Exterior and Interior Regions 92 


a) Far zone radiation from an open-ended 

paral 1 el -plate waveguide: 98 

i) Aperture integration analysis 99 

ii) GTD Analysis 106 

b) Radiation from an open-ended rectangular 116 

waveguide: 

c) Radiation from an open-ended circular 126 

waveguide: 

3.2.2 Transmission of Modal Energy Between 

Two Interior Regions 128 


Transmission of a TEM wave in a paral lei -pi ate 
waveguide into a whispering gallery mode in an 
annular waveguide: 131 


IV 



IV NUMERICAL AND MEASUREMENT RESULTS 132 

a) Reflection from a 2-D horn antenna: 132 

b) Reflection from sharp bends in a 

parallel -plate waveguide: 134 

c) Reflection from a circular bend in a 

parallel -plate waveguide: 137 

d) Electromagnetic backscattering from a 

waveguide cavity model : 139 

i) Results for the Scattering only from 143 

the Rim 

ii) Results for the Interior Cavity Effects 177 

e) Electromagnetic backscattering from a 

circular waveguide cavity: 225 

V CONCLUSIONS 231 

APPENDICES 

A MODAL FIELD EXPRESSIONS IN A RECTANGULAR WAVEGUIDE 235 

B MODAL FIELD EXPRESSIONS IN A PARALLEL PLATE WAVEGUIDE 240 

C MODAL FIELD EXPRESSIONS IN A SECTORAL WAVEGUIDE 243 

D MODAL FIELD EXPRESSIONS IN A 2-D LINEARLY TAPERED WAVEGUIDE 251 

E CIRCUMFERENTIALLY PROPAGATING MODES IN AN ANNULAR REGION 255 

F MODAL FIELD EXPRESSIONS IN A CIRCULAR WAVEGUIDE 265 

G WAVEGUIDE EXCITATION PROBLEM 270 

H EDGE CORRECTION FOR APERTURE INTEGRATION 275 

I APPLICATION OF THE RECIPROCITY THEOREM TO FIND THE D 

RELATION BETWEEN THE SCATTERING MATRICES [S 12 ] AND [S?^ 280 

J EQUIVALENT MAGNETIC LINE DIPOLE 286 

K DESCRIPTION OF THE EQUIVALENCE PRINCIPLE 

EMPLOYED IN THE CALCULATION OF SCATTERING MATRICES 290 

L ILLUSTRATION OF RECIPROCITY IN THE APERTURE INTEGRATION 293 

REFERENCES 297 


v 


LIST OF FIGURES 


FIGURE 



PAGE 

1 . 1 . 

Waveguide geometries considered. 


2 

1 . 2 . 

A waveguide cavity with three junctions. 


7 

1.3. 

Isolation of scattering mechanisms and illustration 
scattering matrices. 

of 

9 

1.4. 

Effective scattering geometry for the problem 
in Figure 1.2. 


10 

2 . 1 . 

An N-port junction. 


14 

2 . 2 . 

An open-ended waveguide cavity problem. 


16 

2.3. 

Scattering by only the rim at the open erd of the 
wavegui de. 


19 

2.4. 

Radiation and reflection at the open end due to a 
waveguide mode incident at the open end. 


20 

2.5. 

Reflection of waveguide modes at junction©. 


20 

3.1. 

The equivalent currents I and M at the rim of the 
open end. ec l ec l 

31 

3.2. 

Diffraction by a wedge. The angles 3 0 , 3 , ip and ip' 
occur in the wedge diffraction coefficient. 

which 

33 

3.3. 

Scattering from the open end of a semi -infinite 
circular waveguide. 


35 

3.4. 

Near field scattering from an open-ended circular 




waveguide using GTD and ECM. GTD, — ECM. 

r=5X, a=lX, e^S 0 , ^-O, <j> r =0 , E 1 ^ 1 


36 

3.5. 

Near field scattering from an open-ended circular 
waveguide using GTD and ECM. GTD, — ECM. 




r=5X, a=lX, e^S 0 , 4 * 1 =0, <(> r =lB0 o , E 1 ^ 1 


vi 



FIGURE 


PAGE 


3.6. Near field scattering from an open-ended circular 

waveguide using GTD and ECM. GTD, — ECM. 

r=5X, a=2X, 9 1 =45°, ^=0, <(> r =0, E 1 =4> 1 

3.7. Near field scattering from an open-ended circular 

waveguide using GTO and ECM. GTD, — ECM. 

r=5X, a=2X, 9 1 =45°, ^=0, <j> r =180°, E 1 

3.8. Near field scattering from an open-ended circular 

waveguide using GTD and ECM. GTD, xxx ECM. 

r=10X, a=5X, 9 1 =45°, (fr^O, /=0, E 1 =<j> 1 

3.9. Near field scattering from an open-ended circular 

waveguide using GTD and ECM. GTD, --- ECM. 

r=10X, a=5X, 9 1 =45°, .MO, ^=180°, t 1 -J 1 

3.10. Near field scattering from an open-ended circular 

waveguide using GTD and ECM. GTD, — ECM. 

r=10X, a=5X, 9 1 =45°, MO, <M90°, E 1 -J 1 

3.11. Near field scattering from an open-ended circular 

waveguide using GTD and ECM. GTD, — ECM. 

r=10X s a=5X, 9 1 =4B°, MO. /=0 ; F 1 =Q 1 

3.12. Near field scattering from an open-ended circular 

waveguide using GTD and ECM. GTD, — ECM. 

r=10X, a=5X, 9 1 =45°, ^O, ^=180°, V 

3.13. Near field scattering from an open-ended circular 

waveguide using GTD and ECM. GTD, — ECM. 

r=10X, a®5X, o' =45°, .J^O, /=90°, E 1 ^ 1 


38 


39 


40 


41 


42 


/ii 


44 


45 


3.14. Scattering from the open end of a semi -infinite 


rectangular waveguide. 46 

3.15. Modal rays in a parallel -pi ate waveguide. 51 

3.16. Modal rays in a rectangular waveguide. 51 

3.17. Modal rays in a circular waveguide. 52 


vii 



FIGURE 


PAGE 


3.18. Modal rays in a 2-D linearly tapered waveguide. 52 

3.19. Modal rays in a 2-D annular waveguide. 53 

3.20. A waveguide junction J between two sections. 54 

3.21. Diffraction of modal rays at the junction discontinuity. 57 

3.22. Excitation of modes by equivalent sources M^ and M d . 59 

3.23. Open-ended parallel -plate waveguide. 61 

3.24. Modal reflection coefficients due to an incident TEjq 

mode in an open-ended rectangular waveguide (a= 2 b). 66 

3.25. Comparisons of the calculated and measured reflection 

coefficients in an X-band and KU-band open-ended 
rectangular waveguide. 67 

3.26. Junction between a rectangular waveguide and a circular 

bend of rectangular cross section. 68 

3.27. A discontinuity in the radius of curvature. 69 

3.28. Reflection coefficient for the TE^ mode in a rectangular 

waveguide due to an E-plane circular bend. 71 

3.29. Reflection coefficient for the TE^ mode in a rectangular 

waveguide due to an E-plane circular bend. 72 

3.30. Junction between two linearly tapered waveguides. 73 

3.31. Junction between two sectoral waveguides. 75 

3.32. Open-ended semi -infinite circular waveguide. 79 

3.33. Modal reflection coefficients due to an incident TEqi 

mode in an open-ended circular waveguide. 82 

3.34. Modal reflection coefficients due to an incident TMqi 

mode in an open-ended circular waveguide. 83 

3.35. Modal reflection coefficients due to an incident TM^ 

mode in an open-ended circular waveguide. 84 

3.36. Junction between paral lei -pi ate and annular waveguide 

sections. 85 


VI 1 1 



FIGURE 


PAGE 


3.37. Reflection of the whispering gallery mode. 89 

3.38. Geometry associated with calculation. 93 

3.39. Open-ended parallel plate waveguide geometry. 98 

3.40. Variation of modal ray angle with frequency. 102 

3.41. Geometry for the edge diffraction analysis. 106 

3.42. Comparison of far zone modal radiation patterns from an 
open-ended parallel -plate waveguide. Mode index = m=5, 

modal ray angle=10°, GTD, AI, xxx modified PTD. 109 

3.43. Comparison of far zone modal radiation patterns from an 
open-ended parallel -plate waveguide. Mode index = m=5, 

modal ray angle=50°, GTD, AI, xxx modified PTD. 110 

3.44. Comparison of far zone modal radiation patterns from an 
open-ended parallel -pi ate waveguide. Mode index = m=5, 

modal ray angle=80°, GTD, AI, xxx modified PTD. Ill 

3.45. Comparison of far zone modal radiation patterns from an 
open-ended paral lei -plate waveguide. Mode index = m= 6 , 

modal ray angle=80°, GTD, AI, xxx modified PTD. 112 

3.46. Comparison of far zone modal radiation patterns from an 
open-ended paral lei -pi ate waveguide. Mode index = m= 6 , 

modal ray angle=85 n , GTD, AI, xxx modified PTD. 113 

3.47. Open-ended staggered parallel plate waveguide. 114 

3.48. Comparison of far zone modal radiation patterns from 
an open-ended staggered, paral lei -pi ate waveguide. 

Modal index = m=3, modal ray angle = 50°, staggering 

angle = t Q =60 o , GTD, __AI, xxx modified PTD. 117 

3.49. Comparison of far zone modal radiation patterns from 
an open-ended staggered, paral lei -pi ate waveguide. 

Modal index = m=3, modal ray angle = 30°, staggering 

angle = t o =60°, GTD, AI, xxx modified PTD. 118 

3.50. Comparison of far zone modal radiation patterns from 
an open-ended staggered, paral lei -plate waveguide. 

Modal index = m=4, modal ray angle = 30°, staggering 

angle = t =60°, GTD, __AI, xxx modified PTD. 119 

0 


IX 



FIGURE 


PAGE 


3.51. Comparison of far zone modal radiation patterns from 
an open-ended staggered, para! lei -plate waveguide. 

Modal index = m=3, modal ray angle = 30°, staggering 

angle = t Q =45°, GTD, AI, xxx modified PTD. 120 

3.52. Comparison of far zone modal radiation patterns from 
an open-ended staggered, parallel-plate waveguide. 

Modal index = m=4, modal ray angle = 30°, staggering 

angle = t Q =45°, GTD, AI, xxx modified PTD. 121 

3.53. A waveguide junction joining two sections. 129 

3.54. Equivalent sources at the aperture radiate in a uniform 

waveguide. 130 

3.55. Termination of a parallel -plate waveguide with an annular 

waveguide. 131 

4.1. The geometry and VSWR of horn antenna. 133 

4.2. A planar 90° bend in a paral 1 el -pi ate waveguide. 134 

4.3. The magnitude of the reflection from the sharp bend of 

Figure 4.2 as a function of bend angle. 135 

4.4. The magnitude of the reflection from the sharp bend of 

Figure 4.2 as a function of section length L. 136 

4.5. A uniform 90° bend in a paral lei -pi ate waveguide. 137 

4.6. Magnitude of the reflection from a circular bend in a 

paral lei -plate waveguide as a function of radius b. 138 

4.7. The geometry of the cavity model. 139 

4.8. Side and top view of the cavity model. 140 

4.9. Geometry of an open-ended rectangular waveguide cavity. 141 

4.10. Radar cross section pattern at f=8.02 GHz, 

“ i a 

E =<t>, <J»=0 plane. 145 

4.11. Radar cross section pattern at f=9.98 GHz, 

A 

E =<)>, <j>=0 plane. 146 


x 


FIGURE 

4.12. 

4.13. 

4.14. 

4.15. 

4.16. 

4.17. 

4.18. 

H IQ 

1 • A ^ f 

4.20. 

4.21. 

4.22. 

4.23. 


Radar cross section pattern at f=ll . 95 GHz, 

E 1 =4>, <(»=0 plane. 

Radar cross section pattern at f=8.02 GHz, 

E 1 =9, <|>=0 plane. 

Radar cross section pattern at f=9.98 GHz, 

“i A 

E =0, <p=0 plane. 

Radar cross section pattern at f=11.95 GHz, 

” i A 

E -9, <f>=0 plane. 

Radar cross section pattern at f=8.02 GHz, 

E 1 = <t» , 4>=90° plane. 

Radar cross section pattern at f=9.98 GHz, 

E 1 =4>, <j>=90° plane. 

Radar cross section pattern at f=ll . 95 GHz, 

E 1= <j>, 4>=90° plane. 

Rsdsr cross ssct 1 on psttorn dt f = 8 > 0 ? GHz* 

E 1 =9, 4>=90° plane. 

Radar cross section pattern at f=9.98 GHz, 

E 1 =9, <J>=90° plane. 

Radar cross section pattern at f=11.95 GHz, 

E 1 =9, <f)=90° plane. 

Variation of radar cross section with frequency. 

E 1 =4>, <(>=0, 9=0. 

Inverse Fourier transforms (i.e., time domain plots) of 
the results in Figure 4.22. 


PAGE 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

159 

160 


xi 



FIGURE PAGE 

4.24. Variation of radar cross section with frequency. 

4> = 0» 9=15° 161 

4.25. Inverse Fourier transforms (i.e., time domain plots) of 

the results in Figure 4.24. 162 

4.26. Variation of radar cross section with frequency. 

E 1 =4>» <t>=0, 9=30° 163 

4.27. Inverse Fourier transforms (i.e., time domain plots) of 

the results in Figure 4.26. 164 

4.28. Variation of radar cross section with frequency. 

E 1 =<*>, <J>=0, 0=45° 165 

4.29. Inverse Fourier transforms (i.e., time domain plots) of 

the results in Figure 4.28. 166 

4.30. Variation of radar cross section with frequency. 

E 1 = 0, <|>=0, 0=0° 167 

4.31. Inverse Fourier transforms (i.e., time domain plots) of 

the results in Figure 4.30. 168 

4.32. Variation of radar cross section with frequency. 

E’-e, 4 >= 0 , 0=15° I 69 

4.33. Inverse Fourier transforms (i.e., time domain plots) of 

the results in Figure 4.32. 170 

4.34. Variation of radar cross section with frequency. 

E 1 = 9, <j>=0, 9=30° 171 

4.35. Inverse Fourier transforms (i.e., time domain plots) of 

the results in Figure 4.34. 172 

4.36. Variation of radar cross section with frequency. 

E 1 = 9, <f>=0, 9=45° 173 


xi i 



FIGURE 


PAGE 


4.37. Inverse Fourier transforms (i.e., time domain plots) of 

the results in Figure 4.36. 174 

4.38. Variation of radar cross section with frequency. 

?-♦. 4>=0. 9=45° 175 

4.39. Inverse Fourier transforms (i.e., time domain plots) of 

the results in Figure 4.38. 176 

4.40. Importance of modes whose modal ray angles are near the 

angle of incidence. 

"i * 

f=10 GHz, E = 4 > , 4>=0 plane. 179 

4.41. Radar cross section pattern at f=8.00 GHz, 

E 1 = i , <{>=0 plane. 181 

4.42. Calculated RCS pattern corresponding to Figure 4.41. 182 

4.43. Radar cross section pattern at f=10.00 GHz, 

E 1 =<f> , <J>=0 plane. 183 

4.44. Calculated RCS pattern corresponding to Figure 4.43. 184 

4.45. Radar cross section pattern at f=12.00 GHz, 

E 1 =<|>, <t>=0 plane. 185 

4.46. Calculated RCS pattern corresponding to Figure 4.45. 186 

4.47. Radar cross section pattern at f=8.02 GHz, 

E 1 =e, <(>=0 plane. 187 

4.48. Calculated RCS pattern corresponding to Figure 4.47. 188 

4.49. Radar cross section pattern at f=9.98 GHz, 

E 1 =9 , $=0 plane. 189 

4.50. Calculated RCS pattern corresponding to Figure 4.49. 190 


xiii 



FIGURE 


PAGE 


4.51. Radar cross section pattern at f=11.95 GHz, 

”i a 

E =0, <|>=0 plane. 

4.52. Calculated RCS pattern corresponding to Figure 4.51. 

4.53. Radar cross section pattern at f=8.02 GHz, 

E 1= 0, <j>=90° plane. 

4.54. Calculated RCS pattern corresponding to Figure 4.53. 

4.55. Radar cross section pattern at f=10.01 GHz, 

E 1= 0, <|>=90 o plane. 

4.56. Calculated RCS pattern corresponding to Figure 4.55. 

4.57. Radar cross section pattern at f=11.96 GHz, 

E 1 =0, <f>=90° plane. 

4.58. Calculated RCS pattern corresponding to Figure 4.57. 

4.59. Radar cross section pattern at f=8.02 GHz, 

E 1= <j>, <f>=90° plane. 

4.60. Calculated RCS pattern corresponding to Figure 4.59. 

4.61. Radar cross section pattern at f=10 . 01 GHz, 

E 1= <(>, <(>=90 o plane. 

4.62. Calculated RCS pattern corresponding to Figure 4.61. 

4.63. Radar cross section pattern at f=11.96 GHz, 

E 1 =<f>, <j>=90° plane. 

4.64. Calculated RCS pattern corresponding to Figure 4.63. 

4.65. Radar cross section pattern at f=10.00 GHz, 

E 1 =4>, <j>=45° plane. 


191 

192 

193 

194 

195 

196 

197 

198 

199 

200 

201 

202 

203 

204 

205 


xiv 



FIGURE PAGE 

4.66. Radar cross section pattern at f=10.00 GHz, 

* 

E =0, <t>®45° plane. 206 

4.67. Radar cross section pattern at f=10.00 GHz, 

E 1 =9, <f>=45° plane. 207 

-■j A 

4.68. Variation of RCS with frequency E =$, <j>=0, 0=0. 208 

4.69. Inverse Fourier transforms (i.e., time domain plots) of 

the results in Figure 4.68. 209 

*• i * 

4.70. Variation of RCS with frequency E =$, <J>=0, 0=15°. 210 

4.71. Inverse Fourier transforms (i.e., time domain plots) of 

the results in Figure 4.70. 211 

4.72. Variation of RCS with frequency E 1 =4», <|>=0 , 0=30°. 212 

4.73. Inverse Fourier transforms (i.e., time domain plots) of 

the results in Figure 4.72. 213 

4.74. Variation of RCS with frequency E 1 ^, <f>=0, 0=45°. 214 

4.75. Inverse Fourier transforms (i.e., time domain plots) of 

the results in Figure 4.74. 215 

-■j A 

4.76. Variation of RCS with frequency E =Q ; d>=0 ; 9=0. o\£, 

4.77. Inverse Fourier transforms (i.e., time domain plots) of 

the results in Figure 4.76. 217 

4.78. Variation of RCS with frequency E 1= 9, <|>=0, 9=15°. 218 

4.79. Inverse Fourier transforms (i.e., time domain plots) of 

the results in Figure 4.78. 219 

“i * 

4.80. Variation of RCS with frequency E =9, <t»=0, 0=30°. 220 

4.81. Inverse Fourier transforms (i.e., time domain plots) of 

the results in Figure 4.80. 221 

4.82. Variation of RCS with frequency E 1 = 9, <J>=0 , 9=45°. 222 

4.83. Inverse Fourier transforms (i.e., time domain plots) of 

the results in Figure 4.82. 223 


xv 



FIGURE 


PAGE 


4.84. The separate returns of modes for the case 

in Figure 4.75. 224 

4.85. Circular waveguide cavity terminated by a short circuit. 225 

4.86. Backscattered field from a circular waveguide cavity, as a 
function of the angle from the axis. All propagating 

modes are included. E =<|> . 227 

4.87. Backscattered field from a circular waveguide cavity, as 
a function of the angle from the axis. Only modes with 
modal ray angles closer than 10 degrees are included. 

E 1= } • 228 

4.88. Backscattered field from a circular waveguide cavity, as a 

function of the angle from the axis. All propagating 

i A 

modes are included. E =0 . 229 

4.89. Backscattered field from a circular waveguide cavity, as 
a function of the angle from the axis. Only modes with 
modal ray angles closer than 10 degrees are included. 

E 1= 0 . 230 

A. l. Rectangular waveguide geometry. 235 

B. l. Parallel -pi ate waveguide geometry. 240 

B. 2. Ray picture of n^ mode. 242 

C. l. Geometry of a sectoral waveguide. 243 

C.2. Projection of ray picture into y-z plane. 250 

C. 3. Projection of ray picture into x-z plane. 250 

D. l. Geometry of a linearly tapered waveguide. 252 

D.2. Ray picture of n th mode. 253 

xvi 


FIGURE 


PAGE 


E.l. Geometry of the annular region. 

E.2. Contours used in the Watson transformation. 

E.3. Location of eigenvalues in the complex v-plane. 

E.4. Ray picture of regular modes. 

E. 5. Ray picture of whispering gallery modes. 

F. l. The circular waveguide geometry. 

F. 2. Conversion of waveguide modal field into the conical ray 

field. 

G. l. Waveguide Geometry 

G. 2. Waveguide geometry with fields radiated to the left 

(n>n + ) of S being of concern. 

H. l. Scattering from a perfectly-conducting half -plane. 

H. 2. The integration path C . 

a 

I. l The geometry of the problem. 

A 

J. l. A wedge illuminated by a S Q -polarized plane wave. 

J. 2. Magnetic line dipoles along z-axis. 

K. l. Radiation of a dipole dp in the presence of the 

semi -infinite waveguide. 

K.2. The equivalent problem of Figure K.l. 

K. 3. The dipole dp g radiating in the infinite waveguide. 

L. l. Modal radiation from the open end. 

L.2. Coupling of incident dipole field. 


255 

259 

262 

263 

264 
266 

269 

270 

273 

275 

278 

280 

286 

288 

290 

291 

292 

293 
295 


xv ii 



CHAPTER I 


INTRODUCTION 


In this report a combination of high frequency and modal ray 
techniques is employed to develop a relatively simple analysis of 
electromagnetic (EM) scattering by electrically large, open-ended 
waveguide cavities with discontinuities. It is assumed that the medium 
exterior and interior to these waveguide cavities is homogeneous, 
isotropic and lossless. The waveguides possess perfectly-conducting 
walls. Therefore, they do not radiate energy transversally, so they are 
in the class of closed waveguides as opposed to open waveguides whose 
cross-section does not confine the field by impenetrable walls. 
Specifically, the two dimensional (2-D) and three dimensional (3-D) 
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different uniformly tapered or curved sections as shown in Figure 1.1. 
Junctions are formed at the open-end and within interior regions where 
different sections are joined. Although the semi -infinite waveguide 
geometries in Figure 1.1 do not indicate any interior termination, the 
present analysis can also account for some simple terminations for which 
the modal reflection matrix is known, such as a planar short circuit, an 
impedance surface termination, or even a homogeneous dielectric loading. 
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a) 2-D Open-ended, flanged, 
staggered parallel plate 
waveguide. 
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b) 3-D Open-ended 

rectangular waveguide 



c) 3-D Open-ended circular 
waveguide. 
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d) Discontinuities formed by 2-D 
linearly tapered waveguide 
sections. 
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Figure 1.1. Waveguide geometries considered. 
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If the waveguide geometry conforms to constant coordinate surfaces of a 
coordinate system in which the Helmholtz equation is separable, then the 
modes can be determined analytically. These analytical modal 
expressions inherently include the variation of the field along a 
non-uniform (or tapered) waveguide as long as modes do not experience 
cut-off as they propagate within the guide. In the waveguide geometries 
considered, the above conditions are assumed to prevail. Therefore, 
under the above conditions, any scattering which occurs within the 
waveguide arises from the abrupt junction discontinuities between 
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different sections. Such a scattering can be expressed in terms of 
modal reflection and transmission coefficients. 

In closed waveguide structures, the modal set is discrete, and the 
classical approach for the analysis of a waveguide discontinuity is to 
expand the fields in terms of the modal sets on each side of the 
discontinuity and then match these fields at the discontinuity by 
enforcing the appropriate boundary conditions. The presence of a 
discontinuity gives rise to a storage of reactive energy in the vicinity 
of the discontinuity in terms of evanescent modes; i.e., the modes which 
decay exponentially with distance away from the junction. Therefore, 
the mode matching approach results in an infinite set of equations to be 
solved. In some cases, exact solutions can be obtained by direct matrix 
inversion techniques, by transform techniques (Wiener-Hopf ) or, by the 
residue calculus technique. In the latter residue calculus technique, 
it is necessary to construct a complex valued function which, when 
integrated along a specified contour gives the infinite set of mode 
matching equations. Also, variational and perturbational techniques are 
used very extensively for waveguides with a single propagating mode. 

The above methods are described in [1,2,3], These techniques are 
limited to certain geometries and they are difficult, if not impossible, 
to generalize to other structures; however, they are usually easier to 
apply to waveguides of smaller cross-sections which allow only a few, 
and in most cases only the dominant, mode or modes to propagate. 

Similar techniques can be applied to open waveguide structures where the 
problem is more complicated by the fact that the modal spectrum has a 
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continuous part. The discussion of continuous transitions in open 
waveguide structures is presented in [4], which also includes a very 
complete list of references on the subject. 

A powerful method to analyze wave propagation through waveguide 
transitions, couplers and tapers is based on the concept of generalized 
transmission line or coupled wave equations. This concept was first 
proposed and studied in detail by Schell kunoff [5]. An extensive 
bibliography and study of different applications of this method are in 
[6]. However, the coupled wave equations are mostly applied to the 
problems of smooth transitions where one converts the Maxwell's 
equations into telegrapher's equations and obtains a matrix differential 
equation for the unknown voltages. A limiting approach can be employed 
for abrupt discontinuities, and in this case, it reduces to the mode 
matching procedure [6]. Adiabatic or slowly varying coupling 
constitutes an extension of coupled wave theory [7], and reduces to the 
WKB approximation. Curved waveguides were studied in detail by Lewin 
[8] and asymptotic approximations with respect to the radius of 
curvature of the waveguide are used for the fields. These approximate 
fields are then matched at the junction to get an integral equation for 
the junction fields. However, coupled wave equations and integral 
equation formulations are related to mode matching procedures; 
consequently, they are cumbersome to employ in practice. 

Application of ray optics to waveguide discontinuity problems is 
relatively new; yet it has received substantial attention in the 
literature because of its usefulness and conceptual simplicity. Both 
Rudduck [9] and Yee, Felsen and Keller [10] applied the geometrical 
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theory of diffraction (GTD) [11] to determine the modal reflection 
coefficients of an open-ended parallel plate waveguide. In [10], the 
effect of multiple interactions between the edges is taken into account 
in a simple, elegant but approximate manner. More rigorous approaches 
to determine the effect of multiple ray interactions at the aperture of 
an open-ended semi -infinite parallel plate geometry were described later 
by Boersma and Lee [12-16], Their results [12-16] agree with the 
asymptotic expansion of the exact solution. An application of the GTD 
to calculate the reflection from the aperture of an E-plane sectoral 
horn was described by Jull [17]. The ray-optical techniques require 
that at each uniform waveguide section, the fields incident at the 
junction be expressed in terms of a set of rays associated with each of 
the characteristic waveguide modes. The latter set is referred to as 
modal rays. For parallel plate, circular and tapered waveguide 
sections, the modal ray representations have been discussed by Felsen 
and Maurer [18,19]. 

Yee and Felsen also described the ray-optical procedure for the 
reflection of acoustic waves from an open-ended circular pipe [20]. A 
ray picture for the EM waves in a circular waveguide is described by 
Ivanyan [21]; however, his rays do not appear to exhibit the usual plane 
wave propagation constant. In [22], a different ray picture which 
consists of a set of conical wavefronts that converge and diverge with 
respect to the waveguide axis is used for the circular waveguide. The 
results based on this modal ray description which are presented in [22] 
are excellent. 


6 


The basic approach employed in this work combines the use of modal 
rays discussed above, with the high frequency techniques, to calculate 
the relevant reflection and transmission properties associated with 
waveguide discontinuities using the multiple scattering matrix (MSM) 
formulation. The general idea behind this approach is briefly 
illustrated as follows. Figure 1.2 shows a waveguide cavity with three 
sections. The scattered field in the exterior region is composed of two 
main contributions; one of these is due to the field scattered from the 
open end alone, and the other is due to the field which is coupled into 
the interior and then re-radiated from the open end. The latter 
undergoes multiple reflections between the open end and the 
discontinuities inside the cavity, after it is initially coupled into 



Figure 1.2. A waveguide cavity with three junctions. 
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the waveguide from the incident field. The fields resulting from these 
multiple interactions may be expressed as a convergent Neumann series as 
done by Pace and Mittra [23]. The same result is directly obtained by 
an alternative procedure based on a self-consistent multiple scattering 
matrix (MSM) method [22]. Therefore, it is possible to identify and 
isolate the wave scattering mechanisms in terms of the scattering 
matrices [S] of the junction which describe the reflection and 
transmission properties of the discontinuity at the junction. In Figure 
1.3 the junction discontinuities are isolated and the wave interaction 
described by the various scattering matrices are illustrated. Once the 
elements of these scattering matrices are determined, then they can be 
combined using the self-consistent MSM procedure indicated above. 
Therefore, the effects of all multiple interactions can be viewed as 
being associated with a single effective discontinuity with its 
corresponding scattering matrix. In particular, the effects of 
scattering from the discontinuities beyond Junction ©for the problem in 
Figure 1.2, can be represented by an effective reflection matrix [S r ], 
as shown in Figure 1.4. 

A dielectric loading inside the waveguide cavity will cause modal 
reflections determined by the geometry and the electrical properties of 
the dielectric material. These reflections can similarly be represented 
by the scattering matrix [S r ], whose entries relate the incident and 
reflected modal field amplitudes. Therefore, Figure 1.4 may also 
simulate the scattering by a dielectric loading located at the 
position ©in waveguide A, as shown in the figure. 
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© 


Figure 1.4. Effective scattering geometry for the problem 
in Figure 1.2. 

As indicated earlier, relatively simple expressions have been 
developed in this work for the elements of the scattering matrices via 
the use of high frequency techniques combined with the modal ray 
description for the fields inside the waveguide. In particular, the 
high-frequency techniques which are employed here include the GTD, the 
equivalent current method (ECM) [24,25,26], and a modification of the 
physical theory of diffraction (PTD) [26,27], depending on the 
scattering matrix being calculated. The previous ray optical approaches 
[9,10] seem to be applicable to certain geometries with some symmetry 
properties. In contrast, the present approach besides being relatively 
simple is more general than the previous ray optical treatments and also 
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allows one to treat many junction types which would otherwise be 
intractable or cumbersome to analyze via the classical mode matching 
approach. Furthermore, since the present approach deals with the 
interactions in a sequential manner, by isolating the effects of the 
different junctions, it is therefore also physically appealing. Some of 
the special junctions treated here have rigorous Wiener-Hopf solutions, 
e.g., open-ended circular pipe and parallel plate geometries [28], 

These Wiener-Hopf solutions are relatively more difficult to obtain and 
are limited to very special geometries such as those mentioned above. 
However, they provide a good check for the present high frequency based 
solutions. As it is common to all ray methods, the procedure described 
here is valid for waveguides which are sufficiently far from cut-off. 
Furthermore, the approximation in the high frequency approaches improves 
as the number of propagating inodes inside the waveguide increases. This 
latter property is actually a merit, rather than a drawback of the high 
frequency based procedures since for small waveguides the other 
approximate methods can safely be used. 

The format of the report is as follows. Definition of the 
relevant scattering matrices and a description of the multiple 
scattering matrix approach which employs these scattering matrices are 
included in Chapter II. Chapter III is devoted to the development of 
scattering matrices for several specific waveguide geometries. They are 
presented under two categories; namely, reflection type scattering 
matrices and transmission type scattering matrices. For a junction 
between two interior waveguide sections, the elements of the scattering 
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matrices are the modal reflection and transmission coefficients. In the 
case of a waveguide discontinuity formed by the open-end; the elements 
of the scattering matrices describe the interior modal reflection as 
well as the exterior rim scattering associated with the open-end; 
likewise, they also describe the exterior radiation and interior 
coupling phenomena. The radiation from the open-ended waveguide into 
the exterior region and the coupling into the interior region due to 
excitation from the exterior may be viewed as being described by 
transmission coefficients. All of these mechanisms are described 
separately, and in some case illustrative examples and comparisons are 
provided. In Chapter IV, the scattering matrices of Chapter III are 
used to solve some two and three dimensional problems. Whenever 
possible, the results are compared with other methods of solution or 
measurements. Finally, some conclusions based on this work are given in 
Chapter V. An ej“t time convention for the fields has been assumed and 
suppressed throughout the analysis. 
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CHAPTER II 


SELF CONSISTENT MULTIPLE SCATTERING MATRIX FORMULATION 

As mentioned in the previous chapter, the solution to the problem 
of the scattering from a waveguide cavity can be constructed by first 
calculating the isolated scattering matrices associated with the 
pertinent discontinuities and then combining these matrices using a self 
consistent MSM formulation. Essentially, MSM is the same as the 
generalized scattering matrix technique ( GSMT ) introduced previously to 
generalize the Wiener-Hopf solutions [2,23], 

The concept of the scattering matrix is widely used in microwave 
circuit theory and relates the scattered wave amplitudes to the incident 
wave amplitudes. Therefore, it characterizes circuit properties of a 
microwave network. 

An N-port junction is shown in Figure 2.1. The incident wave 
amplitudes propagating toward the junction are represented by V + , 

n=l ,N. and scattered wave amplitudes are shown by V“, n=l ,N. 

The scattering matrix [S] of the junction relate these quantities. 

Thus, in general, one can write 
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Figure 2.1. An N-port junction. 


A waveguide containing N "propagating" modes can be represented by 
an N-port junction. The incident and reflected wave (or mode) 
amplitudes are denoted by [V + ] and [V - ]. In addition, one can extend 
the scattering matrix definition to include the scattering of 
"evanescent" modes; however in this case, the modes are not normalized 
to carry unit power. A waveguide cavity with junctions © and© is shown 
in Figure 2.2. It is assumed that the Helmholtz equation is separable 
in the orthogonal curvilinear coordinate system (n,y,5) shown in the 
figure, in which the n-coordinate coincides with the propagation axis of 
the waveguide. Therefore, the n-coordinate corresponds to the direction 
of propagation of the modal field sets inside the cavity. Junction ©is 
the open end of the cavity and junction ©is the termination inside the 
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cavity which is separated from junction 1 by the distance (L), along the 
n coordinate. Therefore, the phase delay of the modal field sets that 
propagate from junction 1 to junction 2 is given by the propagation 
constant of the modes multiplied the distance (L). For convenience, 
only the fields exterior to the waveguide are expressed in spherical 
coordi nates. 

Let a plane wave E 1 be incident at the open end of the cavity. One 
may express incident field as follows: 


"i *i 

E 1 = 0 

n ♦ j 1 e ; 

(2.2a) 

with 



4 - A e 

e-jk" 1 *r 

(2.2b) 


e-jk 1 ' -7 

(2.2c) 

where 

k 1 = -k( sine 1 cos* 1 x + sin6 1 sin<() 1 y + cose^) 

(2.3) 

and 

A 

r = xx 

A A 

+ yy + zz 

(2.4) 


It is noted that k=2ir/A refers to the free space wavenumber, 
(X=wavelength in free space) and r refers to the position vector of an 
observation point. The angles e 1 and are the elevation and azimuth 
angles of the direction of incoming field in the spherical coordinates. 
The corresponding magnetic field H 1 is obtained from 

H 1 = Y Q k 1 x E 1 (2.5) 

where Y = Z is the admittance of free-space. 
oo 
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JUNCTION 

© 


JUNCTION 


Figure 2.2. An open-ended waveguide cavity problem. 


The scattering matrix [Sn] (see Figure 2.3), defines a scattering 
mechanism which converts the incident exterior region field into a field 
scattered back into the exterior region without including any interior 
region effects. In other words, it relates the electric field (E so ) 
scattered from only the open end, to the field E 1 which is incident at 
the open end as follows: 


E S0 


A 1 

X 


E SO 

y 

= [S u ] 

0 

A 1 

E SO 

z 


♦ 





( 2 . 6 ) 


where 


E S0 (P) = 


e*° x + e^° ; 


+ E 


so 


(2.7) 
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and 
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X 0 
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Z 0 
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z<f> 


( 2 . 8 ) 


When the observation point recedes to infinity, i.e., the far zone 
case, then the range dependent part of the scattered field can be 
brought out of the scattering matrix, so [S^] for the latter case is 
defined as follows: 


e!° 


a! 

0 

E so 1 

- [s ii ] 
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A i 

<f> 
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<)> 



with 
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S 9<(. 

s. . 

<p<p— ' 


(2.9) 


( 2 . 10 ) 


The scattering matrix [S 12 ] converts the waveguide modal fields 

incident at the open end (from the interior region) to the exterior 

region fields radiated by these modes as shown in Figure 2.4. The modal 

electric field E* within the uniform waveguide region may be represented 
w 

by 


-+ + - + J ^n 11 

E~ = z C" (e . ± e ) e n 
w n n v nt nn 


( 2 . 11 ) 


Here, e t denotes the transverse (to n) modal electric field of a 

a 

waveguide uniform in the n-di recti on, and & n is the propagation constant 

A 

of the corresponding mode. The axial (or n) component of the modal 
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field is denoted by e nr) . It is noted that the index n denotes a compact 

mode index, which corresponds to the double index nm in 3-D waveguides. 

The superscripts + and - in (2.11) refer to modes propagating in 

+o and -n directions, respectively. It is convenient to define the 

magnetic field H 1 in the waveguide region following the representation 
w 

for the electric field E* in (2.11); thus, the magnetic field is 

w 

given by [1] 


H* = E C* (±h 
w n n L 


. + h ) 
nt nn ; 


* j 5 n n 


( 2 . 12 ) 


The C* in (2.11) and (2.12) are the modal coefficients. If E ro (P) 
denotes the electric field at P exterior to the waveguide region, which 
is radiated by the modes that impinge on the open end, then the 
scattering matrix relates E to E w as follows: 


E ro 

X 


~K°f 

E ro 


[s;°i 
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E r° 

_ Z _ 


J s z°L 



(2.13) 


where the rim is located at n=0, 

E ro (p) = E x ro ; + e™ ; + E ;°; 

and 


CSi 2 ] = 



(2.14) 


(2.15) 
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with 


(2.16) 


t S u°] " t s ui S U 2 S u3 1 

where u corresponds to any one of the cartesian coordinates x, y or z. 
As explained before, for observation points in the far field, the 
scattering matrix [S^] can be written as follows: 



[sj] [C n + ] 



(2.17) 


where 



(2.18) 


OBSERVATION RIM 



Figure 2.3. Scattering by only the rim at the open end of the 
waveguide. 
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Figure 2.4. Radiation and reflection at the open end due to a 
waveguide mode incident at the open end. 
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Figure 2.5. Reflection of waveguide modes at junctionCD. 
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The scattering matrix [S 21 ] describes the transformation or the 
coupling of the incident plane wave field into the waveguide modes as 
illustrated in Figure 2.3. It is clear that [S 21 ] can therefore relate 
to E 1 by 


Cc n '] - cs 21 ] 



(2.19) 


with 


[S n ] - 


tv 

JV. 


( 2 . 20 ) 


The present development can also deal with a near zone source P 

e 

type illumination, in which case the scattering matrix is defined via 
the following relationship: 






where 

ts 2 ;i - [£S X ] [s y ] [S 2 1] 


and 




(2.?A) 


( 2 . 22 ) 


(2.23) 


Here, u corresponds to any one of the cartesian coordinates. 
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It is noted that the problem of determining [S^] 1s re 1 ated to the 
problem of determining [S^] via reciprocity # The precise relationship 
between [S^] and [S^] is d '' scussecl ln Appendix I. 

The scattering matrix [S^] is a modal reflection coefficient 
matrix which is associated with the interaction in Figure 2.4. In 
particular, the elements of [S 22 ] describe the reflection coefficients 
associated with the modes reflected back from the open end into the 
waveguide region when a mode is incident on that open end from within 
the waveguide. Thus, the matrix CS^ 2 ^ relates E“ to E* as follows: 

[C~] = CS 22 1 [C+1 . (2.24) 

The scattering matrix [S r 1, like [S 22 ], is also a reflection 
coefficient type matrix which is associated with the discontinuity at 
junction©. Thus, one can write 

[C+] = [S r ] CC“1 . (2.25) 

As mentioned before, [S ] may also represent a reflection due to a 
dielectric loading, or an effective modal reflection at junction © which 
includes the totality of any other reflection effects arising from the 
presence of all of the discontinuities (or junctions) to the right of 
junction ©in Figure 2.2. 
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At any given operating frequency, the waveguide region can support 
a finite number of propagating modes and an infinite number of 
evanescent (non-propagating) modes. Therefore, the matrices [S^]* 

[S^ ], [$22^ and * CS_] are °f infinite order in a formal sense, however 
one needs to retain only a finite number of elements of these scattering 
matrices in practice, because the distance L shown in Figure 2.2 is 
generally large enough so that the infinite number of evanescent modes 
generated at junction© do not contribute significantly at junction©, 
and vice versa. The finite number of elements of the scattering matrix 
which are retained in practice thus correspond to only the finite number 
of all the propagating (or non-evanescent) modes which can exist within 
the waveguide region. If the distance L in Figure 2.2 is small enough 
so that the "lower order" evanescent modes become important, then one 
must include these "lower order" modes. In either case, in practice, 
one always retains a finite number of elements in the scattering 
matrices which formally may be of infinite order. 

The scattering matrices defined above are combined in a 
self-consistent procedure to obtain the field E s scattered from the 
open-ended waveguide cavity of Figure 2.2 as follows. For an incident 
plane wave E 1 as defined in (2.2), E s can be written as 



3 CR lll 



(2.26) 
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where 


" E x * + + E z z 


(2.27) 


and [R-q] is the effective scattering matrix of the cavity. The 
scattered field E s consists of the field reflected from the rim of the 
open end, and the field transmitted into the waveguide through the open 
end from the incident field at the open end. This transmitted field 
undergoes multiple interactions between the junctions CD and © before 
being radiated from the open end. Let us represent the modal field in 
the waveguide region after taking all these multiple interactions into 
account as follows: 


, + j 8_n 

E' = E T (e . ± 5 ) e 

g m m rnt mn 


(2.28) 


Here, the modal sets are compactly represented by a single index "m" as 
before, and [A + ] and [A"] denote the amplitudes of the modes propagating 
in +n and -n directions, respectively. 

As described earlier, one can decompose the scattered field E s of 


(2.27) as: 



[ s ul 



+ [S 12 ] [A + ] 


(2.29) 


The vector [A + ] represents the amplitudes of the modal field set E* 
propagating in +n direction due to the reflection of the modal field set 
i“ with amplitudes [A“] at junction©. Therefore, one can write: 
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CA + ] * CP] [S r ] CP] CA'] 


(2.30) 


where [P] is the diagonal matrix accounting for the propagation phase 
delay of the modes over the path length (L). 

At junction©, the modal field represented by the vector [A”] is 
excited by the incident field and the reflection of the modes 

A 

propagating in the +n direction from the interior of the open end. This 
effect can be represented by the equation: 


[A“] = [S 22 ] [A + ] + [S 21 ] 


Combining (2.30) and (2.31), one obtains 


[ I_s 22 ps r p -* = [S 2l ] 


Using (2.32) in (2.30) gives: 


CA + ] = [P] [s r l [P] Ci-s 22 PS r P ] _1 [s 21 ] 


(2.31) 


(2.32) 


(2.33) 


Substituting (2.33) into (2.29) yields: 



(2.34) 
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The matrix in curly brackets is the effective scattering matrix [R^] as 
defined in (2.26), so 

[R n ] = [S n ] + [S 12 PS r P [I-S 22 PS/]’ 1 S 21 ] . (2.35) 

As can be observed from the analysis, [R^] includes the effects of the 
presence of the second junction. By this procedure, the self consistent 
method can be directly extended to the analysis of waveguide cavities 
having more than two junctions. 
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CHAPTER III 


DEVELOPMENT OF THE ELEMENTS OF THE SCATTERING MATRICES 
IN THE MULTIPLE SCATTERING METHOD 


In this chapter, expressions for the elements of the scattering 
matrices associated with the discontinuities of waveguide cavities shown 
in Figure 1.1 are developed using a combination of high frequency and 
modal ray techniques. It is noted that there are basically two 
mechanisms described by the scattering matrices: 

1) The reflection mechanism which describes the field reflected 
back into a region due to a field incident on a discontinuity from 
within the same region, and 


2) the transmission mechanism which specifies the tielu coupled 
into a region through a discontinuity when a field is incident on that 
discontinuity from the adjacent region; i.e., from the other side of 
that discontinuity. The reflection and transmission type scattering 
matrices will be treated separately. 


3.1 The Reflection Type Scattering Matrix 

The scattering matrices defining a reflection phenomenon is 
described by either [S^], [$22-1 or [S r ] type scattering matrices. The 
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first one describes the field scattered back into the exterior region 
due to a field incident from the exterior region. The last two describe 
the modal reflection back into the interior region when a waveguide 
field is incident from the same interior region onto a waveguide 
discontinuity. The exterior region scattering (reflection) described by 
the scattering matrix [S^] is going to be discussed separately from the 
case involving interior or modal reflection even though the approach 
employed for both cases are very similar. 

3.1.1 Reflection Back Into the Exterior Region 

In this section, the direct scattering from the rim edge at the 
open end expressed by the scattering matrix [S^] is calculated using 
techniques based on the GTD. According to the GTD, the scattered field 
is initiated from some distinct points (diffraction points) on the rim 
edge as well as from the corners of the rim (if they exist) as a result 
of Keller's generalization of Fermat's principle [11]. In addition to 
the single edge and corner diffracted fields, there exist multiply 
diffracted fields which are produced by rays that undergo multiple 
diffractions across the aperture. These multiple interactions may 
become important if the aperture dimensions are not sufficiently large 
in terms of the wavelength. 

In 2-0 problems the open end is formed by the edges at which the 
semi -infinite waveguide walls are terminated. Therefore, the scattering 
matrix [S^] describes the diffraction of the incident field by these 
edges. The treatment of multiple diffractions is complicated in the 
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situations where one of the edges is in the transition region of the 
shadow boundary of the other edge. j n unstaggered parallel -plate 
waveguide, the edges are exactly on the reflection shadow boundaries of 
multiply diffracted fields. In this special case, the diffracted field 
can be simply decomposed into its ray optical components so that further 
diffractions can be calculated using the GTD [10]. The results of this 
procedure agree with the asymptotic expansion of the exact result up to 
the second order interaction contributions; for the higher order 
interactions, this ray-optical procedure underestimates the asymptotic 
result [29]. However, for numerical calculations the discrepancy is 
negligible. A formally asymptotic ray-optical analysis is proposed in 
[14] based on the uniform asymptotic theory of diffraction (UAT) [30]. 
This lengthy analysis consists of decomposing the interaction field into 
a Taylor series and applying uniform asymptotic theory to each term in 
the series and summing up the results. Another attempt to get an 
accurate multiple diffraction field cunlribuliun is described in 
[12,13]. In this approach, the canonical problem of the diffraction by 
two staggered parallel half-planes is solved and the asymptotic 
approximation to that solution is identified as the modified diffraction 
coefficient (MDC), which includes the presence of the other edge. Later 
on, this solution is extended to non-staggered parallel -plates [15], 

In 3-D problems, the diffraction points migrate around the rim edge 
as the observation point changes position. In some cases, there may be 

a continuum of diffraction points contributing to the scattered field 
which produces a caustic of diffracted rays. The GTD predicts singular 
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fields at a caustic; therefore in order to get a bounded result, one can 
use the so-called equivalent current method (ECM) based on the GTD. The 
ECM is valid provided the diffracted fields are to be found away from 
the shadow boundaries. Away from the caustic regions where GTD is 
valid, the ECM generally blends into the GTD solution as long as the 
waveguide opening is sufficiently large in terms of the wavelength. The 
use of ECM also, automatically, but in an approximate sense, takes into 
account the presence of the corners at the waveguide opening. 

In the ECM, the equivalent currents (I and M ) of the electric 

eq eq 

and magnetic type, respectively, are located at the rim, and they 
radiate in free space to give the diffracted field as shown in Figure 
3.1. The strengths of the equivalent currents are calculated indirectly 
from the GTD, but since they are incorporated in an integral, they give 
bounded results in the caustic regions of the GTD. This definition of 
the equivalent currents makes them correct for the observation points 
lying on the Keller edge diffracted ray cone [11]. In order to extend 
the definition for observation points which are not restricted on the 
Keller cone of edge diffracted rays, it has been proposed to modify the 
expressions heuri stical ly in a manner consistent with reciprocity [31], 
In the following expressions this modification is included. 

The strengths of equivalent electric and magnetic currents are 
given by [26] 
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(3.1) 


and 
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(3.2) 





where l' is the unit vector along the edge direction, and Y , k, E 1 and 
and H 1 are as defined earlier in (2.2) through (2.5). 

/v rs* 

The angles (4», V 3 Q and 3) are associated with the soft and hard 
diffraction coefficients (D $ and D^) which are present in (3.1) and 
(3. 2). In particular, these edge diffraction coefficients 

/V /V 

('I'.V 0 O ,8»<*) in (3.1) and (3.2) are defined for each point on the 
rim, and they are associated with a wedge which locally represents the 


X 



RIM (OPEN END) 


Figure 3.1. The equivalent currents I and M at the rim of the 

. eq eq 

open end. 
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rim geometry. Assuming that Q is any such point of diffraction on an 

fv /V 

edge in an arbitrary curved wedge, the meaning of and 3 with 

respect to the point of diffraction Q becomes clear from Figures 3.2(a) 
and (b) for this wedge configuration. 

The soft and the hard diffraction coefficients are given by [11] 


_i*/ 4 w 

e J sin — 


D. .3 0 » 3 »«) = , / — — 

pi a /2itk /sin 3 sing 


1 


cos a - cos a 


ifl—lp* + 




cos ^ - cos a 


(3.3) 


where the parameter a is related to the exterior wedge angle as shown in 
Figure 3.2(b). 

Finally, the equivalent currents (1^ and M ) are valid provided 
one is observing these sources in directions away from the geometrical 
optics incident and reflection shadow boundaries as mentioned earlier. 
These conditions are certainly met when one is interested in the 
scattering from the waveguide cavity at and near the axial direction. 

The radiated electric field is then given by 

jkZ e -jkR 

E ./ [R X R x I (*') + Y R x M (*')] R-d Z' 

rim e( l 0 e M 

(3.4) 

■ A 

where R=RR is the vector pointing toward the observation point from a 
source point on the rim. 
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(a) (b) 


Figure 3.2. Diffraction by a wedge. The angles 3 0 » 8, ^ and V which 
occur in the wedge diffraction coefficient. 


The use of Equations (3.1) and (3.2) in (3.4) allows one to easily 
identify [S^]. Thus, one obtains 
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and 
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with 



(3.7a;b) 


Examples ; 

a) Near field scattering by the open end of a semi -inf i nite 
circular waveguide: 

The geometry is shown in Figure 3.3. An incident plane wave is 
propagating parallel to the x-z plane (4> 1= 0) and makes an angle (9 1 ) 

A 

with the z-axis. The electric field is evaluated at the observation 
point (P), which is radially separated from the origin by the distance 
(r) and has azimuth and elevation angles <t> r and 0, respectively. The 
results are plotted as a function of 9 in different azimuth planes 
(indicated by <|> ) for different values of r and the radius (a). In the 
near field of the rim of the circular waveguide, there can be four 
points of diffraction for a certain set of observation points (P) as a 
consequence of Keller's generalization of Fermat's principle in the GTD; 
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Figure 3.3. Scattering from the open end of a semi -infinite circular 
waveguide. 

whereas, at other observation points (P) only two points of edge 
diffraction on the circular rim contribute to the scattered field. A 
direct application of the GTD fails in the range of observation points, 
within which such a coalescence of the four points into two points of 
diffraction occurs on the rim edge, since the stationary phase 
condition is violated as these points coalesce. Typically, such a 
disappearance of rays results when three of the four points merge 
together. Conversely, one can also experience a transition from two to 
four rays diffracted to P. A detailed analysis of the migration of 
diffraction points along the rim is given in [32], In such cases where 
the GTD fails, the ECM still gives bounded and reasonably accurate 
results. In the Figures 3.4 through 3.13, calculations based on the GTD 
and the ECM are compared with each other. The comparison gets better in 
the regions where the GTD is valid for large cylinder sizes; elsewhere, 
the ECM is far more accurate. 
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So. 10B 20. 30. 40. 50. 60. 

' 1 THETA (DEGREES) 

Figure 3.4. Near field scattering from an open-ended circular waveguide 
using GTD and ECM. ___ GTO, — ECM. 
r=5X, a=lx, 0^=45°, 4) 1 =0 , <f> r =0, V 
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TOTfiL FIELD (DB) (NORM. TO PIxflxfl) 

-10. -30. -20. -10. 



Figure 3.5. Near field scattering from an open-ended circular waveguide 

using GTD and ECM. GTD, — ECM. 

r=5X, a=lx, o' =45°, if^O, <Ml80°, E 1 
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So. 10. 20. 30. 10. 50. 60. 

’ THETA (DEGREES) 

Figure 3.6. Near field scattering from an open-ended circular waveguide 

using GTD and ECM. GTD, — ECM. 

r=5X, a=2X, e i =45°, ^O, <j> r =0, E 1 *^ 
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• i I I ■ I I 1 1 1 1 I 1 1 1 ' I II* 

So. 10. 20. 30. 40. 50. 60. 

' THETA (DEGREES) 

Figure 3.7. Near field scattering from an open-ended circular waveguide 
using GTD and ECM. __ GTD, — ECM. 
r=5X, a=2X, e 1 ^ 0 , ^=0, <|> r =180 o , f H 
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1 THETR (DEGREES) 

Figure 3.8. Near field scattering from an open-ended circular waveguide 

using GTD and ECM. GTD, xxx ECM. 

r=10X, a=5X, 4. 1 =0, 4> r =0, E 1 ^ 1 
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TOTAL FIELD (DB) (NORM. TO PI*A*A) 

- 10 . - 30 . - 20 . - 10 . 



Figure 3.9. Near field scattering from an open-ended circular waveguide 

using GTD and ECM. _ GTD, — ECM. 

r=10X, a=5X, e 1 ^ 0 , t’-O, ^=180°, E 1 -J 1 
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TOTAL FIELD (DB) (NORM. TO PI*A*A) 

-UO. -30. -20. -10. 



Figure 3.10. Near field scattering from an open-ended circular 

waveguide using GTD and ECM. GTD, — ECM. 

r-10X, a=5\, 9 =45°, * =0, * =90°, E H 
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TOTAL FIELD (DB) (NORM. TO Plxflxfi) 

-IJO. -30. -20. -10. 



Figure 3.11. Near field scattering from an open-ended circular 

waveguide using GTD and ECM. GTD, — ECM. 

r=10X, a=5X, 0 1 =45°, <M0, /=(), E^e 1 
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TOTAL FIELD (OB) (NORM. TO 

-10. -30. -20. -10. 



Figure 3.12. Near field scattering from an open-ended circular 

waveguide using GTD and ECM. GTD, — ECM. 

r=10X, a=5X, e 1 =45°, ^ =0, <t> r =180° , E^e 1 
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10 . 


20 . 30 . 40 . 50 . 60 . 

THETA (DEGREES) 


Figure 3.13. Near field scattering from an open-ended circular 

waveguide using GTD and ECM. GTD, — ECM. 

r=10X, a=5X, e 1 =45°, ^=0, ^=90°, E^e 1 
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b) Near field scattering by the open end 

of a semi -infinite rectangular waveguide: 



Figure 3.14. Scattering from the open end of a semi -infinite 
rectangular waveguide. 


The rectangular waveguide has perfectly conducting infinitely thin 
walls, therefore the diffraction coefficients of Equation (3.3) can be 
written as follows (a=2) 


(♦»♦' »& o ,0) = 


>ir/4 | sin(i|>/2)sin(V/2)l 
2 e“^ \-cos(ip/2)cos(ip'/2)j 1 

cos + * cos *' \/s(nB 0 s1n6 


(3. ft) 


For the incident plane wave fields of (2.2) and (2.5) the equivalent 
currents can be written explicitly as 
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(3.10) 


The equivalent currents, 1^ and M are located on the rectangular 
rim with the dimensions "a" and "b" as shown in Figure 3.14. 

/V 

The angles e Q , e, t* and t in (3.9) and (3.10) are determined as 
ful 1 uwS ; 
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(3.12) 
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The angles t and t are measured from the inner faces of the waveguide. 


cost = - 


( 3 . 13 ) 


cost ~ 


( 3 . 14 ) 


The scattered field at point P with spherical coordinates (r,0,t) due to 
these equivalent currents is given by 
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The result in (3.15) can be simplified if r>>a,b, in which case the 
integrals can be approximated by Fresnel integrals. The numerical 
results for the scattering from the open end of a rectangular waveguide 
are presented in [33] and compared with the Fresnel approximations. 
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3.1.2 Modal Reflection from an Interior Discontinuity 


The scattering matrix which describes the modal reflections due to 
a junction inside the waveguide is examined here using ray optical 
techniques. 

In the ray optical procedure, the incident mode is decomposed into 
ray-optical components, and the GTD is used to determine the field 
diffracted by the junction. These diffracted fields then excite the 
reflected modes. The computation of reflected mode amplitudes has been 
discussed in [10] for the paral lei -plate waveguide geometry using image 
theory. However, this procedure is difficult to apply to other 
geometries which lack analogous symmetry properties. In order to 
overcome this difficulty, it is proposed here to make use of an 
equivalent current approach to get the coefficients of reflected modes, 
which, in the case of paral lei -plate waveguide gives the same result as 
in [10], The proposed procedure can be easily generalized to other 
interior waveguide discontinuities as long as the diffraction 
coefficients for the discontinuities are known. 

Basically, the procedure requires that the modal field expressions 
inside the waveguide sections can be given a ray interpretation. Such a 
ray interpretation of the modal fields for some waveguide geometries in 
terms of the so called "modal rays" is shown in the Appendices. Some of 
these modal rays are illustrated here in Figures 3.15 through 3.19. For 
example, each mode in a paral lei -plate waveguide corresponds to two 
plane waves which propagate in the modal ray directions as shown in 
Figure 3.15. Likewise in the case of the rectangular waveguide, each 
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Figure 3.15. Modal rays in a parallel -pi ate waveguide. 



Figure 3.16. Modal rays in a rectangular waveguide. 
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DIVERGING RAY CONES 



CONVERGING RAY CONES 


Figure 3.17. Modal rays in a circular waveguide. 


CIRCULAR MODAL 
RAY CAUSTIC 


/ 



Figure 3.18. Modal rays in a 2-D linearly tapered waveguide. 
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WHISPERING GALLERY 
MODES 



Figure 3.19. Modal rays in a 2-D annular waveguide. 

mode is decomposed into four plane waves (or modal rays) propagating in 
the directions: 

A ,mr a mir. a 3nm a 

r nm = 1 Cet) x ± (icF) y + nr z • (3.30) 

In a circular waveguide the modal ray picture that is employed in this 
work is shown in Figure 3.17. Each mode consists of two conical modal 
rays which are converging and diverging at the axis of the waveguide. 
Linearly tapered waveguide modes in the 2-D case can be interpreted in 
terms of modal rays bouncing from the two walls. The ray trajectories 
are tangent to the circular modal ray caustic whose radius is determined 
by the mode index as shown in Figure 3.18. It is noted that modal ray 


53 




representation is valid only in the region outside the circular modal 
ray caustic. The 3-D sectoral waveguide modal rays are similar to the 
2-D linearly tapered modal rays, but due to the existence of the third 
dimension the ray trajectories are also oblique to the parallel walls 
(see Figures C.2 and C.3 in Appendix C). Lastly, the ray interpretation 
of 2-D annular waveguide modes is desribed in Appendix E, and they 
include trajectories bouncing from both walls, as well as from only the 
outer wall associated with the whispering gallery modes, as depicted in 
Figure 3.19. 

In order to describe the ray optical procedure which makes use of 
the above modal rays, let us assume that a waveguide section represented 
by A is connected to another section called B, through a junction J, as 
illustrated in Figure 3.20. For convenience, an orthogonal curvilinear 
coordinate system (n,y,S) is assumed in which the Helmholtz equation 
is separable. Furthermore, let the coordinate n coincide with the 



Figure 3.20. A waveguide junction J between two sections. 
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propagation axis of the waveguide section A; it is also assumed that the 
modal field expressions are orthogonal over the waveguide cross-section 
described by the coordinates transverse to n. 

The modal fields in region A propagating through the junction will 
undergo diffraction at the junction and give rise to a set of reflected 
modes in region A which then propagate away from the junction. The 
relationship between the reflected and incident mode amplitudes is 
expressed by the reflection matrix [S^]. The fields entering back into 
the section A as reflected modes are produced entirely via edge 
diffraction of the incident modal rays. Let the n th modal field 


incident at J be represented by 
-+ - - 

E n * < E nt + V * < 3 - 3i 

where, as before, a compact index "n" is used for double modal indices 


and E . is the transverse component of the incident electric modal field 

- ^n^ 

vector; whereas, E is the axial component. Also, e denotes the 


i AA rtU ^ r A f Tlr% A 

y\ upu^juu iuii pnuoc i au uu i • i tic 


propagating in the (±n) direction. 


Let the ray-optical part of the incident modal field undergoing 

diffraction at the junction be represented by E* ... (Note that E + 

r op,m op,n 

is the part of the incident field which is polarized perpendicular to 
the ray trajectory). As explained before, there might be more than one 


modal ray "incident" at a point of discontinuity on the junction for a 
given mode; therefore, the index "i" refers to any one of the "incident" 
modal rays associated with the n t * 1 mode. For example, in the 
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parallel -pi ate waveguide only one of the two modal ray plane waves is 
"incident" at the upper or lower edges, therefore, i=l in that case. On 
the other hand, the four plane wave ray fields in the rectangular 
waveguide are also made up of a pair of rays that propagate toward each 
of the four edges as well as a pair of rays that propagate away from 
those edges. Thus, out of the four modal rays, there exists only one 
pair of modal rays which become "incident" at a given edge; therefore, 
the index "i" goes from i=l to i=2. In the case of a circular 
waveguide, only the axially divergent modal rays are incident at the rim 
edge, so the index "i" takes the value 1. The modal rays in the 2-0 
linearly tapered, and in the 3-D sectoral waveguide are similar to 
parallel -pi ate and rectangular waveguide modal rays, respectively; 
therefore, i=l for the linearly tapered waveguides, and "i" goes from 
i=l to i=2 for sectoral waveguides. Finally, the situation in an 
annular waveguide is clear from Figure 3.19, where the index i=l is used 
in the case of modal rays that bounce from both walls because there is 
only one modal ray that is incident on the discontinuities at inner and 
outer shells, respectively; however, for the whispering gallery modes, 
there is no modal ray incident at a discontinuity on the inner shell 
(and i=0 in the inner shell) and there is only one "incident" modal ray 
for a discontinuity on the outer shell (and thus i=l for the outer 
shell). The corresponding "incident" ray optical magnetic field is 
calculated from: 



A A L 

where r ,j is the unit vector in the direction of propagation of the i tn 
"incident" modal ray associated with the n th modal field. 

As in the calculation of [S^], the ECM will be used to determine 
the entries (or elements) of the reflection matrix [S..]. According to 
the GTD, the diffracted field at a point P, which as shown in Figure 
3.21 is produced by the diffraction at a point Q located at the 
discontinuity is given by 


E d (P) = Z E + . (Q).[e . BD .,B,a(Q)) 

' . . nn.m vv/ ^nm M s 
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| /V 


i=l op,m 


orn s v T m ’ T ’ om 


♦ni * D h (* n i>^ oni »M(Q))] / s(s+p) e " jkS ( 3 - 33 ) 


where a total of N "incident" modal rays are assumed to exist for the 

n th mode. The diffraction coefficients D s and the angles of incidence 

h C: ~ " 

and diffraction are defined as before. The unit vectors 8 . , 8, . 

om m 



CO 


Figure 3.21. Diffraction of modal rays at the junction discontinuity. 
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and t are in their respective increasing angle directions. The distance 
from Q to P is shown by s, and p is the caustic distance for the 
diffracted rays as described in [26]. 

The diffracted field of (3.33) may be viewed as being generated by 
equivalent currents located along the rim of the discontinuity if p»s. 
However, the equivalent currents in this case will be assumed to radiate 
in the presence of a waveguide obtained by the geometric extension of 
the walls of section A, as shown in Figure 3.22. It is noted that the 
equivalent electric edge currents would be shorted out if they were to 
radiate in the presence of the waveguide walls. Therefore, in the 
present work the equivalent edge sources are only of the magnetic type, 
and they radiate within the waveguide. This is different from the 
previous works which employ equivalent electric and magnetic line 
sources in free space [10,14]. 

The equivalent currents are equivalent magnetic line sources and 
magnetic line dipoles given by [34] 


N . * N i /s7_ . D h(*ni , * ;e oni ,p, °^) 

V Q) = l Mj(Q) = -fz VjF (H 0 p n1 -f) 7= 

* i =1 * i=l 0 op ’ ni ? /sing . sing 


(3.34) 
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(3.35) 
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A A 

where %' is the unit vector along the rim of the junction and n is the 
unit vector normal to the walls of the uniform waveguide of Figure 3.22 
at the junction Q) and pointing into the waveguide region. The 
derivation of the expression for equivalent magnetic line dipole is 
Included in Appendix J. The factor 1/2 is included in the expressions 
since the equivalent sources are radiating in the presence of waveguide 
walls. 


<f> 




Figure 3.22. Excitation of modes by equivalent sources M„ and M.. 

x. a 
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The excitation of the modal fields due to these sources can be 

obtained by employing the usual reciprocity theorem as described in 

Appendix G. However, this method must be employed here with a slight 

modification because of the use of "equivalent" currents M. and M. which 

excite the modes; that modification is consistent with reciprocity and 

t h 

so requires that the excitation of the q mode with modal ray angles 
ipqi (i=l,...,N) is due only to the rays diffracted in those directions. 
Therefore, in the expressions for the equivalent sources of (3.34) and 
(3.35) appropriate angles are used to obtain: 
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(3.37) 


Applying the result in Appendix G, one then obtains the expression 
for the reflection coefficients as follows: 


f Hq • (M A +M d ) d*' (3.38) 


^ 2 // ^ qt x V * \ h 5 dYd5 n<m q 
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where h^ and h ? are the metric coefficients (scale factors) for the y 
and 5 coordinates, respectively, and S is the cross sectional area of 
section A. 

The reflection coefficients above give the qn^* 1 element of the 
scattering matrix [S^]. 

The result in (3.37) can be improved for small guide widths by 
including the effects of multiply diffracted rays at the aperture of the 
junction; that effect can be incorporated as an additional term in the 
above expression [10], 

Examples ; 

a) Modal reflection at the open-end of a paral lei -plate waveguide: 





Figure 3.23. Open-ended paral lei -pi ate waveguide. 
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Following the analysis presented in this section and the modal 
field expressions given in Appendix B, one finds that the reflection 
coefficients relating the amplitude of the q th reflected mode to the 
amplitude of the mth incident mode are explicitly given by 

i ) TM z case; 



where e . $ and 6 for the n th mode are as defined in Appendix A and 
om m n 

the angles and are shown in Figure 3.23. Also shown in the 
figure are the waveguide width (a) and the staggering distance (d). 


ii) TE z case; 



(3.40) 


The results in (3.39) and (3.40) agree with that obtained in [10] via a 
different procedure. 
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b) Modal reflection from the open-end of a rectangular waveguide: 

The geometry and the modal field expressions are presented in 
Appendix A. The elements of the reflection coefficients which relate 
the pqth reflected mode amplitude to nmt* 1 incident mode amplitude are 
calculated via (3.38) and given explicitly as follows: 
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( 3 . 43 ) 

( 3 . 44 ) 

( 3 . 45 ) 

( 3 . 46 ) 

( 3 . 47 ) 

( 3 . 48 ) 

( 3 . 49 ) 

( 3 . 50 ) 

( 3 . 51 ) 

( 3 . 52 ) 
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N , Y , n , and m. are as defined in Appendix A. 
nm’ nm* a* b 

The result in (3.41) is plotted for different incident and 
reflected mode values by Dr. C.D. Chuang of the Ohio State University, 
ElectroScience Laboratory and presented in [351. Some of those results 
are shown in Figures 3.24 and 3.25. In Figure 3.25, the calculation is 
compared with measurements [36], In this case, the waveguide dimensions 
allow only the dominant TEjq mode to propagate. It is noted that, even 
though at such low frequencies the neglected multiple order diffractions 
at the open end become important, the comparison between the measurement 
and calculation is quite reasonable. One would expect the comparison to 
get better as the frequency increases. 
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50 


0 



measured calculated (Equation (3.41)) 

Figure 3.25. Comparisons of the calculated and measured reflection 
coefficients in an X-band and KU-band open-ended 
rectangular waveguide. 
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c) Reflection of TEjq mode from an E-plane circular bend in a 
rectangular waveguide: 

x 

a 

b 

TE (0 MODE 

Figure 3.26. Junction between a rectangular waveguide and a circular 
bend of rectangular cross section. 



The dimensions of the rectangular waveguide shown by "a" and "b" 
allow only TE^g mode to propagate. This rectangular waveguide is joined 
to a uniform circular bend of rectangular cross section, as shown in 
Figure 3.26. The mean radius of the circular bend is shown by R. In 
this case the edge diffraction coefficients in the expressions for 
equivalent sources given by (3.36) and (3.37) will have to be modified, 
because the diffraction is now due to a discontinuity in the radius of 
curvature as shown in Figure 3.26. One can employ the diffraction 
coefficients presented in [37] for this case; they are expressed as 
follows: 


c +i'fr/4 / 2 

Mvj.v,,) = e V« k < F * G > 


(3.58) 
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0 


Figure 3.27. A discontinuity in the radius of curvature, 
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where 


a2-ai l+cos(vi+V2) 
W (cosvi+cosv2) d 


(3.59) 


F = - 


and 


G 


a2-aj l+cos(vj-V2) 
2k (cosvi+cosv2) J 


(3.60) 


In (3.59) and (3.60), and a^ refer to the curvature of the surfaces 
which make up the discontinuity, Vj is the angle of the incidence and 
is the angle of diffraction, measured from the tangent direction 

A 

(y-direction) as shown in Figure 3.27. 

The expression in (3.41) for the reflection coefficient can be 
applied here together with the diffraction coefficients of (3.58) to 
calculate the reflection of the TE 10 mode. This calculation is done and 
plotted in Figures 3.28 and 3.29 (denoted by GTD) as a function of 
radius R for two different values of the waveguide height "b" . This ray 
optical result is found to be independent of the waveguide width "a" and 
is not varying appreciably with the value of "b" for comparably large 
values of R. A result based on the variational approach for the 
equivalent circuit parameters is presented in [38], The approximate 
reflection coefficient based on that approach in [38] is also presented 
in Figures 3.28 and 3.29. That reflection coefficient is likewise very 
insensitive to the dimension "a" for the selected range of parameters. 
However, its magnitude increases rapidly with the dimension "b" which is 
unexpected based on the physics of the problem. 
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i 


Figure 3.29. Reflection coefficient for the TE^q mode in a rectangular 
waveguide due to an E-plane circular bend. 
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d) Modal reflection from a junction between two linearly tapered 
waveguides: 



Figure 3.30. Junction between two linearly tapered waveguides. 


The reflection coefficients for this 2-D configuration can be 
determined employing the modal field expressions of Appendix D in 
(3.38). The amplitude of the reflected mode due to an incident 
mode is given by 


i ) TE X case: 



/M~ e" jir/4 

/ 6 £ 2<t> 

on oq o 



W 9 lq * g o"f g "f g l ) 

/ p i 


- j C B n ( p 2 ) +0q ( P 2 ) ) p 2 + J C ®n ^ P 1 ^ + 6 q ^ P 1 ^ p 3 

P2 


~ 7T ~ 7T 


D h (9 2n» 9 2q ; V2» 6= 2» a 2 ) 

(3.61) 
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ii ) TM x case: 




/27k e 


•j*/ 4 


qn " 2* 


(_!)<!♦" °s( V V 8 o'f s ~f' , l1 
/ BnlPjlS (P^ P! 


-j( 6 p ( P 2 )+0 ( p 2 ) ) p 2 + j( p l^ + 0 q^ p l^ p ] 


/i^yopi) p2 




(3.62) 


In (3.61) and (3.62), the parameters e Qn , 6 n , <|) o and are defined 

in Appendix D. The angles a^, and are determined from the wedge 
angles as shown in Figure 3.30, and the modal ray angles 9 ln and e 2n for 
the n th mode is calculated via the expression given in Figure D.2, for 
edges® and©, respectively. 
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The modal field expressions have been given in Appendix C for TE X 
and TM X type of modes. The ray interpretation of these modes is also 
discussed in the same appendix. The reflection coefficients are 
presented below for TE X and TM X incident modes. 

i) TE X case: 

In order to excite the pqt* 1 reflected mode, equivalent magnetic 
line sources due to an incident nm th mode are located at edges® and © 
and they are given by 


M 


*1 

%2 


. 2 
x Z 
i=l 



L Y n (k2-(H1)2) 
nm O'- a ’ 

<W P 1> 



Kr^xd-i) 


• JS nm (o l )p l + Js nm (p l )p l 
2 2 


D. ( 9i X ,0i X ;3 ,3 ,a. ) 
h v £nm* £pq ’ n p r 

rj 


(3.63) 


The upper characters correspond to edge ©and the lower ones give the 

value of the equivalent source at edge©, as illustrated in Figure 3.31. 

The modal ray angles (corresponding to nm th mode) 0 x nm and 0^ present 

in the diffraction coefficients are determined via the expression given 

in Figure C.2 for edges© and©, respectively. Finally, the acute angle 

nn 

6 is determined using cose = 177 • 
n n 
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The amplitude of the reflected pqth mode due t0 an nm th i ncident 
mode is obtained as follows: 


pq,nm 


nr~ ✓ kMi 1 ) 2 /(k2-(f L )2 

* V 2jk 77 ; 7 7 — A 


^ ^ e on e om e op e oq ^o 


np no' 


(_j)q+m ^h^ 9 lnm* 0 lpq’ 6 n*^p* a l^ 


/0 nm (p l )0 pq (p l ) P 1 


^> 2 >% q <p 2 >VJ(U»i>Vn»»i x x - ( ~ , 

+ h v 2nm* 2pq’ n* p* 2\ 


nm v w 2 ' pq v w 2 ' v 2 


(3.64) 


ii ) TM X case: 


In this case, the equivalent magnetic line dipoles are located at 
edges ® and © and are given by: 


M 


ft,? 

1*“"* a f ~J 


dl 

d2 


c /Sir u nm 

2 2 


|(- 1 H • , , ,,mr, , 

t(-l) i “ J 0 nm VP l' p l' i ' J '- a J*' 1 " 1 ' 
2 e 2 2 


) 

" [p V 1 -(kiT) + X (lea) cose^ pq ] (3.65) 

a 

For the modes away from cut-off, the x component of the magnetic 

a 

dipoles are small compared to the radial (p) component. Therefore, in 
the calculations it is assumed that the equivalent line dipoles have a 
radial orientation and excite a pq th reflected mode due to an incident 


D s (9 lnm» 9 lpq ;0 n ,0 p ,a l 
2 2 


4 sine 


ip« 
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nnr n mode. It is noted that the equivalent currents of (3.65) can 


excite a TM^ mode as well as a TE x mode. Therefore, the reflection 
coefficient corresponding to a pq*"* 1 type TM x mode will be represented by 

pp t h 

R ; whereas the reflection coefficient corresponding to a pq type 
pq,nm 

he 

TE mode will be represented by R„ . The explicit expressions for 
x pq,nm 

pp he 

R and R are presented below: 
pq ,nm pq,nm 


R pq,nm “ “VJF <(> 0 


8tt qir \A 2 “(a ) 2 


np on 


r k 2 -(|p) 2 ^ e on e om e op e oq ^o 


(-l) q+m 


D (9? ,9? ;B ,3 ,a, ) 

s v lnm’ lpq’ n’ p’ l 7 


l/^nmW^pqW) p f sin9 lpq 
j ^nm {p 2 )+B pq (p 2^ p 2 +j ^nm (p l )+J ' 0 pq (p l^ p l 
\/ 8 nm^ p 2 )0 pq^ p 2 ) p 2 sin9 2pq 


D s ^ 9 2nm * 9 2pq ’ 0 n * 0 p ’ a 2 ^ 


(3.66) 
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pq VM r s v lnm’ lpq’ n* p’ 1 

. i\ m+q 

(_1) 2 sin9 lpq ^nm(Pl) 9 pq(Pl) Pi 

‘ j ( 0 nm ( p 2 ^ +0 pq * p 2^ p 2 +j ( 0 nir/ P 1 ^ +0 pq ( P 1 ^ P 1 


0 pq 0 2nm’ 0 2pq ;0 n’ 0 p’ a 2^ 

2 sin9 ^pq V^mMW^ >2 


(3.67) 



f) Reflection from the open end of a circular waveguide: 


The geometry of the circular waveguide is shown in Figure 3.32. 

The modal expressions and their modal ray interpretation are left to 
Appendix F. 

i ) TE z case: 

The equivalent currents and at the rim of the open end 

excite TE z and TM z modes. The reflection coefficient corresponding to a 

reflected TE, (TMl mode will be represented by R®® (r|)® 1. They 

z z pq,nm pq,nm 

are given by 


X 



Figure 3.32. Open-ended semi-infinite circular waveguide. 
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* ee 

pq;nm 


'n'kZ N N J (P* ) J (P * ) 
o pq nm n pq n nm 

j4a cos6 pq + cos6 nm 

2 2 


[n cos <S cos6 »f 

L p pq nm 


k a t nn sin6^ sins «g] 5 nn 
on pq nm 3 np 


(3.68) 


2 N N nJ (P* )J'(P ) 

. irk pq nm n v pq' n v nm' 

R n = TJ : 1 1 [cos6 sin6 »f 

J /»rtC X J_ r/\e i L r»/i r» m 


pq,nm 


COS6 pq + COS6 nm 


pq nm 


+ sin6 cos 6 »gl <$ 
pq nm 3 np 

(3.69) 


where 


^pq ^nm AD 

f = cos — cos — [U2 (sec 5pq ♦ sec« nm ) 


(3.70) 


6 5 

pq nm AC 

9 " sin — sin — [1-2 (secS pq + secS m ) 


(3.71) 


A - (-1)' 


-j(2ka-ir/4) 


4/irka 


(3.72) 


C = 1 + l (jB)*/ A FT 
*=1 


(3.73) 


D = 1 + I (-jB) / /t+T 
*=1 


(3.74) 


B = (.l) n Ye’ j2ka 


(3.75) 
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It is noted that the effect of all the multiply interacting rays 
across the aperture are included in the above expressions for the 
circular waveguide opening; in particular, that information is contained 
in the "f" and "g" terms given above. 


ii) TM z case: 


Similar to the TE„ case, the reflected modes may be TM or TE 
z J z z 

type. The reflection coefficient corresponding to a reflected TM z (TE z ) 

mode will be represented by nm (R®j] ), The expression for R eh 

pq,nm pq,nm' K pq,nm 

is the same as (3.69), namely; 


R eh = R he 
pq,nm pq,nm 


(3.76) 


and 


? hh 

pq ;nm 


wkY N N J'(P )J'(P ) 

o pq nm n v pq 7 n v nm' 2 2 

T7 ; fk a r sin<$ <;inS »f 

j4a cos 6 nn + cos6 nm - on pq nm 


pq 


- n cos 6 cosfi »g] 6 
pq nm 3 np 


(3.77) 


The results in (3.68) through (3.77) are compared with the exact 
Wiener-Hopf solutions for various modal reflections at the open-end of a 
circular waveguide. The comparisons are shown in Figures 3.33 through 
3.35. The ray-optical solution closely approximates the exact solution 
especially at higher frequencies. 
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present solution Wiener-Hopf [28] 

Figure 3.34. Modal reflection coefficients due to an incident TMgi mode 
in an open-ended circular waveguide. 






g) Modal reflection in an annular waveguide terminated by a 


parallel -pi ate waveguide: 


E* 



Figure 3.36. Junction between parallel -pi ate and annular waveguide 
sections. 


In this example, the reflection coefficients corresponding to 
annular waveguide modes are derived by considering the discontinuity at 
the junction between an annular and a paral lei -plate waveguide as shown 
in Figure 3.36. The annular waveguide has outer and inner shell radii 
shown by "a" and "b", respectively. The polar coordinates p and $ are 
used to define a point with respect to the center of the annular 

A 

section. The goemetry is assumed to be infinite in z-direction, and in 
this example only the reflection of TE z modes will be considered. The 
treatment of TM z modes follows similar lines. The modal field 
expressions in the TE z case are obtained from a magnetic field given by 


H* = + z Ajj (ka,kb) f*(kp)e n 
so that 

-+ a 1 3 -+ » 1 3 -+. 

E n " p jkp Y Q H H n “ *** jk Y Q 3 p H n 
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(3.78) 


(3.79) 



In (3.78), the function f^(kp) is given by 


= u(2)' 


f ”(kp) - H^' (kp) 


.( 1 ), 


.(I)' 


H' 7 (kp)-H' 7 (kp) 
n v n 


H (2) (kp) 


(3.80) 


p=b 


p=b 


where the primes denote differentiation with respect to the argument, 

and the v 's are determined from 
n 


df S 


= 0 


p=a 


(3.81) 


h 

The normalization factor (ka, kb) is used in (3.78) such that 
for a propagating mode one obtains 


1 

2 



E + x (H + )* dp = 1 
n ' n 7 


(3.82) 


To determine the normalization constants, one can generate the Green's 
function as done in Appendix E for TM^ polarization. Alternatively, one 
can also use the following identity [39]: 


/ w v (kp)W v (kp) 


dp kp • 3W 3W\ 

p = -2 V l* v 9 V “ w v 3 V J 


(3.83) 


where w v (kp) and W v (kp) are combinations of cylindrical functions (i.e., 
Bessel, Neumann, or Hankel functions) of order v and primes denote 
differentiation with respect to the argument. 
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The equivalent magnetic line sources for the excitation of the q th 
reflected mode are located at E and E 1 of Figure 3.36 and their values 
are given by 



V V 



A* (ka.kb) 
2 


-H^ ' (kb) H^(ka) 
v n v n 

H^'(kb) H^(kb) 


D h^ 5 ln ,5 lq ) 


2 2 
(3.84) 


where D^j is the diffraction coefficient for the discontinuity in the 
radius of curvature (such as the one in (3.58)), and the incident and 
reflected modal ray angles are as determined in (E.31) and (E.34). The 
subscript 1 or 2 is relevant to the equivalent source at edge 0 or @ of 
Figure 3.36, respectively. 

The reflected q** 1 mode is obtained from (3.38) such that 
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For whispering gallery modes, the line source Mg corresponding to 
the discontinuity at edge ©on the inner circle (of radius=b) has to be 
discarded since the whispering gallery modal rays do not illuminate this 
discontinuity. Therefore, the general expression for the reflection 
coefficient in (3.85) reduces to the following for whispering gallery 
(WG) modes: 



(kb)H, (2) (ka)D^(6 1n ,5 


In’ lq ; 


(3.86) 


The result in (3.86) for the magnitude of the reflection of a 
whispering gallery mode is plotted in Figure 3.37 as a function of the 
radius "b" while keeping (a-b)=0.4x (wavelengths). The diffraction 
coefficient of (3.58) is used in the calculation. For the range of 
radius "b" shown in Figure 3.37, it is found that the only propagating 
mode is the whispering gallery mode. The modal ray direction of the 
whispering gallery mode was sufficiently away from grazing (more than 20 
degrees); therefore, the diffraction coefficients of (3.58) could safely 
be used. One observes that the amplitude of the reflection coefficient 
decreases as the radii of inner and outer circles are increased. 
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3.2. The Transmission Type Scattering Matrix 


The scattering matrices [S^l and [S 21 ] characterize the 
transmission of the field through the open end by relating the 
coefficients of modal field inside the cavity to the amplitude of the 
exterior field. One may also encounter the transmission of modal energy 
from an interior region to another adjacent interior region through the 
junction between them. These transmission mechanisms are discussed 
separately in this section. 

In order to obtain expressions for the transmission coefficients 
one may be tempted to use the GTD equivalent currents employed in 
calculating the reflection coefficients. However, it is noted that the 
GTD equivalent current concept fails at and around the geometrical 
optics shadow boundaries where the diffracted field is not ray optical. 
This was not a problem in the calculation of reflection coefficients, 
because the incident and reflected modal ray angles were away from the 
shadow boundaries for the modes not close to cut-off. However, in the 
case of transmission through a junction, the modal rays corresponding to 
the transmitted modes may be close to the shadow boundaries of the 
incident modal rays, and their associated reflected modal rays, at the 
discontinuities. Therefore, to obtain the transmission coefficients, 
one can integrate the incident modal field across the aperture as in the 
Kirchhoff-Huygens 1 approximation for aperture integration (AI). 

Although the asymptotic evaluation of this integral gives some end point 
contributions to account for the edge effects, they are not as accurate 
as predicted by the GTD [40]. Therefore, in this section a modification 
of the PTD presented in [26] will be employed to get a proper correction 
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to the edge effects predicted by the Kirchhoff -Huygens' approximation. 

The PTD approach discussed in [26] modifies the original PTD developed 

by Ufimtsev [27] by employing equivalent currents to refine the 

approximate physical, optics (PO) integral for the fields scattered by a 

conducting body. These equivalent currents are again placed at the 

edges of the scatterer, and their expressions are similar to the GTD 

equivalent currents as in (3.1) and (3.2); however, the GTD diffraction 

coefficients D s are now replaced by Ufimtsev diffraction coefficients 
h 

D<j which are given by; 

D r v f 

PO 

In (3.87), D <5 is the PO diffraction coefficient obtained from the 
asymptotic end-point contribution of the PO integral. Therefore, the 

PTD, in this format, requires an integration of the geometrical optics 
(GO) currents over the surface of the scatterer (which is the PO 
approximation) and another integration of the Ufimtsev type equivalent 
currents over the edges of the scatterer. The latter tends to correct 
for the incomplete edge effects contained in the PO approximation. The 
important property of the Ufimtsev type equivalent currents is that, 
unlike the GTD equivalent currents, they can be employed for observation 
points in the transition regions of the shadow boundaries. It is noted 
that when the result of PO integration is identically zero, then 
Ufimtsev diffraction coefficients have to be replaced by the GTD 
diffraction coefficients. One thing still remains to be shown; namely, 
the same Ufimtsev type equivalent currents employed to correct the PO 
approximation for the scattering by a perfectly-conducting surface with 
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an edge also corrects the incomplete edge effects contained in the 
Ki rchhoff-Huygens 1 approximation for the AI considered in this work. 

This is done in Appendix H where it is noted that the Ki rchhoff-Huygens ' 
integral basically plays the same role as the PO integral. 

3.2.1. Transmission Between Exterior and Interior Regions 

The interior to exterior coupling is the result of the radiation of 
modal energy from within the waveguide to the external free space region 
through the open end, and it is described by the scattering matrix 
[S^]. The reciprocal problem of transmission into the waveguide region 
via a coupling of the exterior field incident at the open end is 
characterized by the scattering matrix [S^]. As might be expected, 
[$2^] can be directly related to [S^] via reciprocity; this relation is 
discussed in Appendix I, and [S^] therefore be found here directly 
from [S 2 P]. 

The scattering matrix 3 describes the transmission or coupling 
of the incident plane wave field into the waveguide modes as illustrated 
in Figure 2.3. [S,^] i s defined in (2.19). 

For later convenience, let the waveguide excitation be an electric 
current moment dp g which lies at the point P in the region exterior to 
the waveguide as shown in Figure 3.38. [S 2 ^] will be developed here 

for estimating the coupling of the fields of the exterior source dp g 
into the interior waveguide region, via the open end. This more general 
situation reduces to the special case of plane wave incidence on the 
waveguide opening as in Figure 2.3 if the source dp g is allowed to 
receed to infinity. 
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03 



a) Part of H due to coupling of the direct field of dp g 
through the aperture. 



b) Part of [S^] due to diffraction correction from the rim. 


(17*0 PLANE) 



c) Equivalent problem. 

Figure 3.38. Geometry associated with [S^] calculation. 
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The field E~ which is coupled into the waveguide region by dp in 

W G 

Figure 3.38 may be expressed as in (2.11) by 


E" = Z C" (E . -E ) 
w n nt nn 
n 


+j6„n 


= E C" E" 
^ n n 


(3.88) 


The associated magnetic field H~ is likewise given via (2.12) as 

w 


- - - + J' B n n 

H' = Z C (-H . +H ) e n = SC H 
w n nt nn „ n n 

n n 


(3.89) 


The modal coupling or transmission coefficients (C~) in (3.88) and 
(3.89) may be found from the equivalent problem in Figure 3.38(c) which 
illustrates an equivalent surface and line source distribution at n=0 
within an infinite waveguide which is a geometric extension of the 
semi -infinite waveguide of Figure 3.38(a,b). The equivalent sources in 
Figure 3.38(c) generate the same fields in the waveguide region as those 
in Figures 3.38(a) and (b) if the equivalent sources are found exactly, 
as described in Appendix K. Here, the equivalent sources are determined 
from asymptotic high frequency techniques so that the field coupled into 
the waveguide are approximations to the true fields therein. These 
approximations are high frequency approximations which are expected to 
work well even down to the lowest propagating mode in the waveguide. 
According to the high frequency estimates based on the GTD, one obtains 

3. ~ 3f + 3^ (3.90) 

and 

M s ~ M s G ° + M* . (3.91) 
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" pA -pA 

where j“ and represent the unperturbed or goemetrical optics (GO) 

field produced by dp g within the aperture region but in the absence of 

the waveguide structure. The additional contributions to and M s must 

arise from the diffraction by the edges of the aperture which are 

denoted by and in (3.90) and (3.91). It is easily seen that 5®° 
~T0 

and M 5 may be expressed by 

3f = [n x Hd Pe ] (3.92) 


and 




(3.93) 


at the waveguide opening. The unit vector n defines the unit normal at 

the aperture surface pointing into the waveguide. In (3.92) and (3.93) 

Edp e and Hdp e are the electric and magnetic fields incident at the open 

-i 

end from the external source dp g at P. The incident fields (E<jp e anc ^ 

-i 

Hdp e ) represent the unperturbed fields of dp & which exists in the 
absence of the waveguide. Clearly, the incident fields are given by 



jkZ 0 M 
~TT s 


x s 



(3.94) 




(3.95) 


The expressions in (3.94) and (3.95) are valid for distances s 1 
which correspond to dp g being in the near zone of the waveguide 
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aperture. However, s 1 cannot be made extremely small to where the 
reactive field terms of the type l/(s ) and 1/ ( s ) become significant; 
these higher order range dependent terms are ignored in (3.94) and 
(3.95) which pertain only to the radiation fields. Before proceeding to 
calculate the modal transmission coefficeints (C~), it is convenient to 
decompose (C~) as follows: 


C" = C G0 + C d 
n n n 


(3.96) 


where (C GI ? is the part of (C - ^ due to J GG and M GG ; whereas (C d \ is the 

' nJ ' n> s s ' n> 

part (C~) which is produced by J d and M d . It is now an easy matter 
rn — pn “Pn 

to find (C p ) from j” and M™ by employing the results of Appendix G 
which indicate the manner in which the electric and magnetic current 
sources excite modes inside a waveguide. Thus, from (G-8) of Appendix 
G, the GO part of the transmission coefficients are given by 


.GO 


■ If [i + • j w - H + * m go ] ds 1 
aperture n s " s 


2 / / 

aperture 


§♦ * «r 

n n 


ds 


(3.97) 


At this point it is also worth mentioning that the AI method, which 
is used to calculate the GO part of transmission into the waveguide and 
radiation from the waveguide to the exterior region, satisfies the 

reciprocity principle; this fact is shown in Appendix L. The ed 9 e 
contributions (C d ) may be calculated via the Ufimtsev equivalent edge 
sources given by (for N incident modal rays of each mode as in (3.33)): 


N -i - N I U/E I D h ^ *^ni ;g o ,g nj) (hV) 

M * = ^ L- Y o *J k J sine o sin3 ni 2 


(3.98) 
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and 



M 


« 


. N 
n E 
i=l 



D s ^ *^ni ’ e o ,0 ni ) 
l sin *nil /sin3 0 sine n . 


(E 1 •*' ) 
— 2 — 


(3.99) 


In (3.98) and (3.99), the incident ray from dp g makes an azimuthal angle 

A 

V about the unit tangent i' at any point of diffraction on the rim. 

(see Figures 3.2(a) and (b)). Likewise, ij» . in (3.98) and (3.99) 
corresponds to an azimuthal angle which the ray diffracted from the edge 

« A 

along the s^ direction into the waveguide region makes about V at the 

A* 

same point of diffraction [note; e m . and are like g and i|> in Figures 
3.2(a) and (b); likewise, g Q and ip 1 also have the same meaning as in 
those figures], A more accurate analysis will include the obliquity of 
the equivalent magnetic dipole source of (3.99) with respect to the 

A 

axial direction (n). However, this effect is negligible for the modes 
sufficiently far from cut-off; hence, it is not shown in (3.99). The 
coefficients (C d ) excited by and are explicitly given by 


C d - 
n 


N 

E 

j=l 


/ 

rim 

2 I / 


ft • ft 


mJ)] 


dt' 


K x 


H + **ds ' 
n 


(3.100) 


The scattering matrix [S 2 j] can easily be identified via (3.96), (3.97) 
and (3.100). In the case of an infinitesimal electric current moment or 
dp e type illumination as indicated in (3.94) and (3.95) the scattering 
matrix can be identified via the following relationship: 


■ [ [ s x„] cs yn ] [S zn ] ] 





(3.101) 
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where 


t S 21 ] * t IV ^ 1 


dp ex 


* 

x' • dp 
K e 

dp 

ey 

= 

A 

y'* dp e 

dp 

_ ex_ 


z‘ • dp 
_ e_ 


(3.102) 


(3.103) 


If the source receeds to infinity, then Ep e * E 1 > corresponding to a 
locally plane wave illumination. In this case [Sgj] takes the form of 
[$2^] given in (2.19) and (2.20), by bringing out the range dependent 
part. 


Example: 

a) Far zone radiation from an open-ended paral lei -plate waveguide: 



Figure 3.39. Open-ended paral lei -plate waveguide geometry. 
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i) Aperture integration analysis 

As shown in Figure 3.39, the geometry is independent of the 

a 

y-coordinate and the width of the waveguide is shown by "a". The AI 
analysis of the TM^ case will be considered first; the corresponding TE y 
case will be discussed briefly later on in a similar fashion. The 
expressions for the surface equivalent currents at the aperture are as 
fol 1 ows : 


J 


eq 


and 


n x H' 



(3.104) 



A 

x n 



(3.105) 


A A 

where n=z is the unit normal vector pointing into the free space at the 
surface of the aperture. 

The far field radiated by these equivalent currents is given by 


(The distance from the origin to the far zone observation point is shown 
by p, in Figure 3.39.) 


£rad 


(p,9) = y 



e " JkP ,/JiT » 

r r dx - y 8 tt / (“ cose ) 

p o 


p-jkp 

E 1 — pr dx 1 
x /p 


(3.106) 

or 


E rad (p,e) 



, a imr 
/ sin X 

o 


8 -J kp 

/p- 


dx' 


cose + 



(3.107) 
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The integral in (3.107) can be evaluated and the result is given by 


j^rad 


(p,e) 



. / .mir 2 

cose +Vl-(ki'J 


-jk[p-fsine] 

e 


2j 



2 sine + 2 ? 


r ka mir, 

sin sine - -?] 

k mir 

2 sin9 ' 2 ? 


(3.108) 


As seen from (3.108), the far field pattern of each mode is 
sinx , 

composed of two ~ ^ — (sine) functions with their peaks in the modal ray 
directions. Further simplification occurs in the expression if one 
separately examines the odd and even modes. For m=2n+l, where n is an 
integer one obtains 


I cose +\A-(fj)' 


E rad = ^C m VS [ cos e + 

,ka . 

2mir cos [~2 sine) 

■ o a A 2 
sin20 -(ka) 


-jk[p-?sine] 


/F 


(3.109) 


and for m=2n, one finds that 


i rad . cj& 


cose + 


V 1 ,mir. 2 
l-(k?) 


-jk[p-2sine] 


/ p * 


,ka . ^ 

2mir sin (“2 sine) 

jFI . r mir 2 
sin20 -(k?) 


(3.110) 
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It is noted in the above equations that if one considers the 

radiation only from odd (or even) modes, their peaks occur at sine^^f 

for the m th mode pattern. In this direction, the other odd (or even) 

ka 

modes have a null in their pattern due to the cos( — jr sine) (or 
kd 

sinf - ? sine)) factor in the numerator. Therefore, around this m th modal 
ray direction, the pattern is mainly governed by the m th mode. 

The aperture integration is reasonably accurate in the front 
half-space (z>0) for frequencies exceeding the cut-off frequency by 5% 
[28]. It is also known that AI gives exact radiation in the modal ray 
angle directions [28]. In Figure 3.40, the variation of the modal ray 
angles with increasing frequency is plotted for several odd modes. 

Modal ray angles start at 90° from the waveguide axis at cut-off and 
sharply decrease with slight increase in frequency. In fact, the slope 
of the curves in Figure 3.40, are unbounded at the cut-off frequencies. 
As the frequency is increased above the cut-off, the modal ray angle 
decreases monotonical ly ; however, the rate of decrease becomes slower. 
Therefore, at a fixed frequency (or ka), the angular separation between 
the different modal ray directions becomes smaller for points near the 
axis of the waveguide than farther away from the axis, as can be 
observed in Figure 3.40. 

For all modes, at a frequency 5% above the cut-off the modal ray 
angles are about 71° from the axis. The accuracy of the aperture 
integration is limited to modes with modal ray angles less than 71°. 

For frequencies very close to cut-off, the AI result becomes unbounded 
which is not only inaccurate but also nonphysical. In [28], it is shown 
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from the exact Wiener-Hopf solution for this problem that the radiation 
for frequencies very close to cut-off is very small. 

As mentioned before, the form of (3.108) is such that the radiation 
in a fixed direction (9) is mainly governed by the contributions of 
modes whose modal ray directions are closest to 9. Therefore, if the 
odd modes with modal ray directions in the 9^ neighborhood of 9 are 
included, then the error resulting from neglecting the contribution of 
the other odd modes can be calculated as follows. 

ka 

Let 9j=9-9^ and 9 2 s9+9 b* Also one may define Nj=c~sin9^ j and 
ka 

* sina^ 3 where the symbol t • 3 denotes the odd integer part for 
the odd modes. From (3.109), the error resulting from the exclusion of 
odd modes with ray angles outside the 9 b neighborhood of 9 is therefore 
given by 


R 


o 


Ni . N . 

E 1 E Pad + E E rad 

n=l n=N ? +l 


(3.111) 


where N=Ctpsin9 j, in which 9. < 90° but 9„ is close to 90°. 
i i a 

Therefore, in (3.111), it is assumed that the contribution of only a few 
(at most a couple of) modes with modal ray angles greater than 9 is 
neglected due to the fact that AI must be limited to modes not too close 
to cut-off for reasons indicated above. 

From Equations (3.109) and (3.111), one finds 



4_ 

< /8irkp 



n=l 

n=odd 


C„ A„ sine, 
n n l 


N 


C n A n sin9 t 


sin 2 9-sin 2 9 + J sin 2 9-sin 2 e o 
1 r^N^+l 2 

n=odd 


(3.112) 
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where A is the amplitude of the n th odd mode. For backscattering 
n 

problems, these modes are excited by an incident plane wave with an 
incident angle. In Appendix I, it is shown that the amplitudes A^ of 
these excited modes can be determined from the radiation problem via 
reciprocity and they are given by 





ka 

~2 sin9 nrr 


/' ka - 

cos (“2 sine) 

sin 2 e - (If) 1 


(3.113) 


Inserting these mode amplitudes into (3.11?) one obtains 


sin 3 e 1 (sine^-sinQg) sin 2 ©^ 

cos0j [sin 2 0-sin 2 0j ] cos 0^ [sin 2 0-si n 2 0 2 ] 

(3.114) 

A similar analysis can be performed for the even modes of (3.110). 
The important thing to be noticed from (3.114) is that, as the frequency 
increases, the number of propagating modes increases also; however, by 
retaining only the modes with modal ray angles in the 0^ neighborhood of 
0, the error remains bounded. Also in the region close to the axis 
(i.e., 0+O) one can reduce the value of 0^ and still expect to get the 
same error. This is due to the fact that modal ray directions are more 
densely clustered around the axis. 


1 

I /?°l < ItTr? 
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For the TE^ case, the equivalent currents are written as follows: 


and 


- A “*i A *i 

j = n x H 1 = -X H 
eq y 


% “ E x n = y E; 


(3.115) 

(3.116) 


The radiated field is given by 

JH 

Y 8ir 


H Pad = J 


-jkp 


a . e “jkp j— a e ' 

-y o VS o / Ex ~7v~ dx ' "V87 / ( “Cos 9) H 1 , 7= dx' 


V /p 


(3.117) 


or 


_ -jkp 


“ 

, a mir e 

J cos a x dx 


\ /. ,mir.2 
cose + VI- (lea) 


Qrad * r v \/IK 

H “ yC m Y «iV8w 


(3.118) 

The integral in (3.118) can be performed in closed form and the 
result is given by 


H rad = y C Y V §7 
J m m OTr 


\!H 

.v a it 


. ) ,mir 2 
cos e +\/ 1- (kF) 


e 


/p* 2 


,-m 


ka mw 

sin [ 2 sine + 2J 


ka mir " 

m sin | 2 sine - ~o\ 

J ~ k — iff + (-j) m 


7 sine + 


k . mir 
7 sine - ■ja 


. (3.119) 


Again, the pattern is composed of two ~ — type functions. An 
error analysis similar to the TM case can be performed in this TE case 

y y 

to indicate the effect of excluding modes whose modal ray angles lie 
outside the e^ neighborhood of 6. 
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ii) GTD Analysis 


In the GTD approach, the far field radiation is due to the 
diffraction from the edges of the half planes making up the waveguide as 
shown in Figure 3.41. This analysis is described below for the TMy case 
using only a first order GTD (which neglects rays multiply diffracted 
across the aperture). The analysis for the TE y case is similar, hence, 
it will not be discussed here. 



Figure 3.41. Geometry for the edge diffraction analysis. 
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(3.120) 


For the TM y case. 


E 1 = C 


m 


mu mu 
e -e 

2j" 


-ja z 
m 


as before. The modal plane wave corresponding to the first term on the 
right hand side of the previous equation is incident on the bottom edge 
at (x=0,z=0) and the other one corresponding to the second term is 
incident on the top edge at (x=a,z=0). 

The electric field which is diffracted from the bottom edge at 
x=0 is given by 


d 1 e- Jkp C m (-l)e-J* /4 

E yl ~ ^y D s /p" /p" 2 /2irk 


1 


ir-9-\p_ 

COS 


•FT— 9+\p_ 

cos C— 2 — 


(3.121) 


using 


e ; 


r 

y 'm 

ij 


and 


rriTT 


sm + n = Yi 


(3. 122a, b) 


t h 

where is the rn modal ray direction (see Figure 3.41). 

Similarly, the diffraction from the top edge at x=a, can be written 


as 


e -jkp e -jkp / i\ e 

i „ c e jmiT e jkasin9 1 L 


-jif/4 


E y2 = E y D s /p“ = /p 


2 /5ttTT 


ir+0— \p 

cos ( m 


ir+9+\p 


■) cos (") 


(3.123) 
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The total diffracted field is obtained by superposing the 


diffractions from each edge; thus, one obtains 


E = 


E d , + E d 

yi y2 


-jkp 




M. 

, V 8ir 


1/ a 

j - 2Sine 


T 


•m 


, ka rmr 

sin [ 2 s i n9 + 2 ! 

0+lp m 

sin (-2-3 


(-J) 


m 


ka rmr " 

sin [^2 sin 9 - “ 2 ] 


9-4 


. (3.124) 


• / an 

sm (_-) 


The result in (3.124) can be improved by including the contribution 
to the radiation resulting from the rays which undergo multiple 
diffraction across the aperture [10], However, for wide apertures 
(approximately greater than one wavelength) these higher order effects 
are small enough to be neglected. 

The Equations (3.108) and (3.124) based on AI and GTD, 
respectively, have the same limiting values at 0=+^ [41]; however, away 
from the modal ray directions, the AI result differs from the GTD 
result, and as described in Appendix H, the AI result can be corrected 
using a modified PTD approach so that it agrees with the GTD. 

In Figures 3.42 through 3.46 the pattern factors of certain modes 
obtained using AI, GTD and modified PTD are compared. As can be seen, 
the modified PTD approach gives results almost indistinguishable from 
the GTD result. Also, the agreement between the GTD and AI gets better 
for larger apertures. 
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Figure 3.42. Comparison of far zone modal radiation patterns from an 
open-ended parallel -plate waveguide. Mode index = m=5, 
modal ray angle=10°, GTD, AI, xxx modified PTD. 
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Figure 3.43. Comparison of far zone modal radiation patterns from an 
open-ended parallel-plate waveguide. Mode index = m=5, 
modal ray angle=50°, GTD, AI, xxx modified PTD. 
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Figure 3.44. Comparison of far zone modal radiation patterns from an 
open-ended paral lei -plate waveguide. Mode index = m=5, 
modal ray angle=80°, GTD, AI, xxx modified PTD. 
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Figure 3.45. Comparison of far zone modal radiation patterns from an 
open-ended paral lei -plate waveguide. Mode index = m=6, 
modal ray angle=80°, GTD, AI, xxx modified PTD. 
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The above results can be extended to the calculation of radiation 
from a pair of staggered, semi -infinite parallel plates as shown in 
Figure 3.47. The staggering is determined by the angle (t 0 ). The 
analysis is developed only for TMy case here as an illustration. In 
this case, the equivalent currents of (3.104) and (3.105) take the 
following form 

5 eq = y (-cost Q H^ + sint^j) (3.125) 

and 

M eq = (-z cost Q + x sint Q ) . (3.126) 



Figure 3.47. Open-ended staggered parallel plate waveguide. 
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The radiated field due to these equivalent sources is given by 


i; ad (p,9) 


- _ .M 
y c m V 8 tt 


a 7 mir 2 cos(t o + 0 )‘ 
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k 
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* ■ • \ w u 1 
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Equation (3.127) reduces to (3.108) for the non-staggered case when 

IT 

t 0 =2 as expected. 

The first order GTD result of (3.124) becomes 


■jkp 


E = —f=— C 
y /p m 




l“ a , 4 v/ J.,2 COS(t o +0)- | 

Jk |_2 cot(t 0 ) V l-(ka) “ sin(t 0 )_| 
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m sin(aA) m sin(aB) 

9+4» " 0-4) 

sin( HI.) sin( HI) 
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127) 


128) 
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Figures 3.48 through 3.52 show the comparisons of the results 
obtained by AI, GTD and the modified PTD approaches. In all of the 
cases, the results based on the modified PTD agree with the 
corresponding GTD results. 

The multiple diffractions between the edges can be treated as in 
[10] and can be added to the first order GTD result. This is done and 
compared with other methods for the transmission problem in [42], 
Therefore, the higher order diffraction analysis is omitted here. 


b) Radiation from an open-ended rectangular waveguide: 


The far zone radiated field for the open-ended rectangular 
waveguide geometry of Figure A.l can be obtained in closed form such 
that 


jjjrad 


5 E 9 + ♦ % 


(3.131) 


where E 0 and E^ can be further decomposed into contributions associated 
with the Kirchhoff-Huygens* (or aperture) integral and the equivalent 
Ufimtsev edge current integral. Thus, the electric field components are 
given by 


E 9 “ E 9k + E 0U 


and 


E. = + E a 

<f> 4>k 4>u 


(3.132) 

(3.133) 
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Figure 3.48. Comparison of far zone modal radiation patterns from 
an open-ended staggered, paral lei -plate waveguide. 
Modal index = m=3, modal ray angle = 50°, staggering 
angle = t =60°, GTD, AI, xxx modified PTD. 
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3 # 49. Comparison of far zone modal radiation patterns from 
an open-ended staggered, parallel-plate waveguide. 
Modal index = m=3, modal ray angle = 30°, staggering 
angle = ^=60°, GTD, AI, xxx modified PTD. 
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3.51. Comparison of 
an open-ended 
Modal index = 
angle = t =45° 


far zone modal radiation patterns from 
staggered, parallel -plate waveguide. 
m=3, modal ray angle = 30°, staggering 
— GTD, AI, xxx modified PTD. 
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The contribution from the Kirchhoff Huygens' approximation is given by 


• o +m ^ ^ j ^ 
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for TE mode 
for TM mode 




for TE mode 
for TM mode 


(3.134) 


(3.135) 


(3.136) 


(3.137) 


(3.138) 


(3.139) 
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Likewise the contributions from the equivalent Ufimtsev edge currents 
are given by 
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0u V 8 tt j nm r 
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and 


123 


where 
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(3.141) 


u, v, A + , B + are defined as before, and D s is given in (3.42), and 


f» nm /k 
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124 



(3.144) 


sine oh = Cl - (n g /k) 2 ] 1/2 

7 7 1/2 

sin0 h = (1 - sin 9cos «|>) (3.145) 
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(3.146) 


(3.147) 

(3.148) 


and 

sin0 v = Cl - sin 2 esin 2 <|>] 1 ^ 2 . (3.149) 

As seen in (3.134), (3.135), (3.140) and (3.141) the pattern has a 

sinx 

x “ type behavior in both the 0 and $ directions. There are four 
sinx 

different ~ forms, each of which is due to the integration of a plane 
wave component of the mode as described in (A. 16). Therefore, each ~~ y X 
exhibits a peak at the corresponding plane wave direction. Also, as the 

waveguide dimensions get larger electrically, the main beams get 
sharper. As a result, only a few modes contribute strongly to the 
radiation in a given direction of observation. 

Numerical results for this case are left to Chapter IV, where the 
backscattering analysis of a waveguide cavity model is performed and 
numerical values are compared with experimental measurements. 
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c) Radiation from an open-ended circular waveguide: 


For the open-ended circular waveguide geometry of Figure F.l, the 
far zone field can be written as 



-sinne 

cosne 


) + J E 


4> 


(Cosn <f>\ -| 
''Sinn<|>J 



(3.150) 


where Eg and E^ can be separated as contributions from the 
Ki rchhof f-Huygens 1 approximation and the equivalent Ufimtsev edge 
currents such that 
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The contribution from the Ki rchhoff -Huygens 1 approximation becomes 
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™nm modes incidence : 

E ek “ J kN nm P nm 2(cos6 nm -cos 6) J n^ P mJ J n ( kasin0 ) (3.156) 


E 



(3.157) 


Likewise, the contribution from the equivalent Ufimtsev edge currents is 
gi ven by 
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^nm modes incidence : 
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where 


f( 0 ,i ) 

nm 


cos 


- cos 


nm 

T 


cos5 nm - cose 


(3.162) 


The validity and accuracy of these results are extensively 
discussed and numerical results are presented in [22], 


3.2.2 Transmission of Modal Energy Between Two Interior Regions 


A waveguide junction between two waveguide sections is shown in 
Figure 3.53. It is of interest to determine the transmission 
coefficient T qn which gives the coefficient of the q th mode transmitted 
into guide "B" when an n th mode is incident on the junction from within 
the guide "A". These transmission coefficients T qn then are the qn th 
entry in the scattering matrix [Sba]. The T qn is found via the modified 
PTD approach involving an aperture integral and a Ufimtsev type 
equivalent current integration as discussed previously in 3.2.1. The 
equivalent currents at the aperture are calculated in terms of the 
incident field in the Ki rchhoff-Huygens ' approximation for the aperture 
integral. In addition, the Ufimtsev type equivalent currents which 
represent a correction to the above aperture integral approximation are 
located at the aperture edges; i.e., at the edges of the junction 
forming the aperture. Furthermore, the surface equivalent currents at 
the aperture are given by 


and 



n x H 


n 




A 

x n 


(3.163) 

(3.164) 
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Figure 3.53. A waveguide junction joining two sections. 


A 

where n is the unit vector at the aperture S. pointing into region "B" 

a 

as shown in Figure 3.53. These equivalent sources radiate into region 
"B". 

An application of the reciprocity theorem to an appropriate set of 
fields in region "B" (as shown in Figure 3.54) will yield the strength 
of the modes transmitted into region "B". The Ufimtsev type equivalent 
currents are determined from the ray optical parts of incident electric 
and magnetic fields represented by n - and respectively, and 

they are given by 


1 a / 57 
V 3F 
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Figure 3.54. Equivalent sources at the aperture radiate in a uniform 
waveguide. 


where the angles \j>qi and 3q-j are computed from the i th ray trajectory of 
the qth mode in region "B". Such an application of the reciprocity 
theorem was indicated earlier in Section 3.1, and hence it will not be 
described here. For the reasons presented earlier, only the n-directed 
component of the magnetic dipole source is included in (3.166). 

The transmission coefficients are then given via reciprocity 
arguments by 
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(3.167) 
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Example: 


Transmission of a TEM wave in a paral lei -plate waveguide 
into a whispering gallery mode in an annular waveguide. 

The geometry is shown in Figure 3.55. The dimensions and 
polarization allow only a TEM mode to propagate in the paral lei -plate 
waveguide, and only one whispering gallery mode can likewise exist in 
the annular waveguide. The coupling from the incident TEM mode into the 
whispering gallery mode is found to be very close to unity for the 
values of radius "b" changing from IX (wavelength) to 2X. 


* TEM MODE 





Cb*0.4X 


WG '^CIRCULAR 
MODE PORTION 


Figure 3.55 Termination of a paral lei -pi ate waveguide with an annular 
waveguide. 
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CHAPTER IV 


NUMERICAL AND MEASUREMENT RESULTS 

In this chapter, the scattering matrices described above are 
combined to solve some two and three dimensional problems of interest. 
Whenever possible, the results are compared with other methods of 
solution or measurements. 

a) Reflection from a 2-D horn antenna: 

The voltage standing wave ratio (VSWR) of the two-dimensional horn 
antenna is calculated. The antenna is shown in Figure 4.1 and has a 
waveguide width of 0.375 wavelengths which is chosen so that only the 
dominant TEM mode will be excited far from the edges. As seen from the 
figure, the reflection is due to the discontinuities at the throat and 
the open end. The problem is solved using the procedure described in 
Chapter III, and the result is compared with the moment method 
calculation [43] in Figure 4.1. In the latter figure, the reflection 
due to the throat alone is shown separately, and as the horn angle (a) 
gets smaller, the reflection from the throat gets smaller also. 

However, for o small, there is strong reflection due to the 
discontinuity at the open end. 
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b) Reflection from sharp bends in a parallel-plate waveguide: 

The reflection of the TEM mode due to sharp bends in the parallel 
plate waveguide of Figure 4.2 is calculated. The reflection and 
transmission coefficients are determined from the formulas developed in 
the previous chapter. The magnitude of the total reflection from both 
junctions is presented in Figures 4.3 and 4.4 for various bend angles 
(9) and lengths (L) of the bent section, respectively. 



-INCIDENT TEM WAVE 
► REFLECTED TEM WAVE 
( FROM BOTH JUNCTIONS) 


Figure 4.2. A planar 90° bend in a parallel -plate waveguide. 
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REFLEC 



Figure 4.3. The magnitude of the reflection from the sharp bend of 
Figure 4.2 as a function of bend angle. 




REFLECTION COEF 



SECTION LENGTH (WAVELENGTH) 


Figure 4.4. The magnitude of the reflection from the sharp bend of 
Figure 4.2 as a function of section length L. 
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c) Reflection from a circular bend in a parallel -pi ate waveguide: 

The geometry of the circular (90°) bend in a parallel -pi ate 
waveguide is shown in Figure 4.5. The magnitude of the total reflection 
from both junctions is calculated as a function of the radius (b) as 
shown in Figure 4.6. Note that the reflection in the present case is 
very small compared to that which occurs in the case of a sharp bend as 
shown previously in example b). 


.<////</////. 


/T1 0.4X ~ 

q *>>>** ' rrr . , ^ 


~ INCIDENT TEM WAVE 
►REFLECTED TEM WAVE 
PROM both junctions) 




V 

V 

s> 
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Figure 4.5. A uniform 90° bend in a parallel -pi ate waveguide. 
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d) Electromagnetic backscattering from a waveguide cavity model: 


PERFECTLY 



Figure 4.7. The geometry of the cavity model. 


The waveguide model of interest in this work is shown in Figure 
4.7. It is basically an open-ended cavity composed of two waveguide 
sections. The first section is part of a sectoral waveguide with one 
end open; whereas the other end of this section is connected to a second 
section which is a uniform waveguide with a planar termination at its 
far end. The exterior of the second section is curved at the back end 
to minimize the scattering coming from the exterior features of the 

A 

structure. The axis of the waveguide coincides with the z-axis of the 
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coordinate system, and the y-axis is in the vertical direction. The 
model is made of wood and then coated with a conductive paint. The 
dimensions of the cavity are shown in Figure 4.8. 

The side walls in each of the waveguide sections are parallel to 
the y-z plane; therefore, there is no tapering effect to be included for 
those walls. The first step in the analysis is finding the modal field 
expressions in the rectangular and sectoral waveguide sections. This is 
done in Appendices A and C where the expressions are given, and the 
relationship between the mode sets of sectoral and rectangular 
waveguides is shown. Therefore, the effect of tapering can be included 



(a) SIDE VIEW 



(b) TOP VIEW 

Figure 4.8. Side and top view of the cavity model. 
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by modifying the propagation constant in the rectangular waveguide mode 

A 

expressions which are transverse to the x-direction. So, the analysis 
is first done on a rectangular cavity as shown in Figure 4.9, and the 
solution is modified subsequently to include the effect of tapering. 

The geometry of the rectangular waveguide structure is shown in 
Figure 4.9. The cross-sectional dimensions of the waveguide are "a" and 
"b" in the x and y-coordinate directions, respectively. The length of 
the waveguide from the open end to the back wall is given by the 
dimension (L). All of the walls of the cavity are assumed to be 
perfectly conducting. The structure is illuminated by an incident plane 
wave which is given by 


:i /pia . c il\ JMxsine^os*)) 1 + ysi n 0 1 cos <j> n + zcose 1 ) 
: = (E 0 e + e 


(4.1) 


where 0 < 9 1 < ir/2, 0 < $i < 2ir are the elevation and aspect angles of 
the incident field direction, respectively. 


y 



Figure 4.9. Geometry of an open-ended rectangular waveguide cavity. 
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The scattered field is composed of two main contributions: 1) the 
field scattered from the open end by itself, and 2) the field which is 
radiated from the open end. The latter undergoes multiple reflections 
between the open end and termination at the back wall after it is 
initially coupled into the waveguide from the incident field. However, 
from experimental measurements, these multiple interactions are 
determined to be negligibly small for the cavity model being considered 
here; therefore, their effect is ignored in this study. Only the first 
order interaction is discussed which includes the coupling of the 
incident energy into the interior waveguide modes through the open end, 
and subsequent reflection of these modes from the back wall and finally, 
the radiation of these reflected modes from the open end. In Appendix 
I, it is shown that the mechanisms of coupling into and radiation from 
the open-end are equivalent via the reciprocity principle. 

For a finite cross-section, the field incident at the open-end 

A 

excites a finite number of modes which propagate in the -z direction 
without attenuation, as well as an infinite number of evanescent modes 

A 

which attenuate exponentially away from the open end in the -z 
direction. In this study, it will be assumed that the dimensions "a" 
and "b" are large enough to excite at least a few propagating modes and 

the length "L" is long enough for the effects of the evanescent modes to 
be negligible. 
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i) Results for the Scattering only from the Rim 


To check the validity of the analysis, the scattering from the rim 
is calculated using (3.11), and results are compared with a set of 
measurements obtained using the Compact Range Facility at the 
ElectroScience Laboratory, The Ohio State University. The measurements 
were performed on a cavity model as shown in Figures 4.7 and 4.8. In 
order to remove the interior cavity effects, the inner surface of the 
back wall of the cavity was covered with absorbing material. Two sets 
of measurement results were obtained. 

a) At a fixed frequency, the angle (9) from the z-axis is 
varied in the horizontal ( 4>=0) (or vertical (<|>=90 0 )) plane from 0 to 90 
degrees. 

b) At a fixed angle (9) in the horizontal (i>=0) (or vertical 

( <J>=90°) ) plane the frequency is varied from 8 GHz to 12 GHz in 10 MHz 
steps. Therefore, a bandlimited frequency response is obtained. From 
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illustrate the scattering mechanisms. This is accomplished by 
processing the frequency response through a Kaiser-Bessel window and 
inverse Fourier transforming the windowed result using an FFT algorithm. 
Since the measured spectrum is bandlimited, the time domain response 
represents the impulse response of the target which is convolved with 

2 s1n V 

7 l cosu) Q t (4.2) 


where co^half -bandwidth (2 GHz), and <o o is the center frequency which is 
10 GHz in this case. 
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The measured and calculated radar cross section (RCS) patterns 
are shown in Figures 4.10 through 4.21 for different polarizations 
and frequencies. In each figure, measurements (solid line) and 
calculations (dashed line) are drawn on the same scale in terms of dB 
relative to square meter (dBSM). Figures 4.10 through 4.15 correspond 
to a horizontal (x-z) scan; whereas. Figures 4.16 through 4.21 
correspond to a vertical (y-z) scan. In both the horizontal and 
vertical scans, the radar cross section is measured and analyzed for two 

A A 

polarizations (<j> and 0). Finally, each case is repeated at three 

different frequencies, namely around B, 10 and 12 GHz. In both the 

* 

horizontal and vertical scans, a directed incident electric field will 
not diffract from two edges which are parallel to the scan plane. 
Therefore, the contribution to the radar cross section in these cases is 
dominated by the diffracted fields from the remaining two edges. As 
shown in Figures 4.10 through 4.12 and 4.16 through 4.18, the 
calculations agree reasonably well with the measurements. However, for 

A 

the other polarization (0), all four edges of the open end contribute to 
the radar cross section, and the contribution from two of these four 
edges reduces essentially to that from the end points (or corners). As 
shown in Figures 4.13 through 4.15 and 4.19 through 4.21, the agreement 
between the measured and calculated results is not as good as the 

A 

^-polarization case, especially in the plane (<f>=0) pattern, horizontal 
The reason for this discrepancy will become clearer after discussing the 
frequency domain responses, suffice it to say for now that it is due to 
the imperfections of the model. 
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RADAR CROSS SECTION (DBSM) 



LT>0. 30. 60. 90. 

' THETfl (DEG) 


Figure 4.10. Radar cross section pattern at f=8.02 GHz, 

E 1 =<}>, (j»=0 plane. 

measured, — calculated 
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90 . 


Figure 4. 


60. 

THETA (DEG) 


. Radar cross section pattern at f=9.98 GHz, 

— j * 

E = 4 > , <j)=0 plane. 

measured, — calculated 
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Figure 



.12. Radar cross section pattern at f=11.95 GHz 

E 1 = 4 >, <t>=0 plane. 

measured, — 


cal culated 


90 . 


60. 

THETA (DEG) 


Figure 4.13. Radar cross section pattern at f=8.02 GHz, 

E 1 =9, 4>=0 plane. 

measured, — calculated 


RflDRR CROSS SECTION (DBSM) 



Figure 4.14. Radar cross section pattern at f=9.98 GHz, 

E 1 =0, f=0 plane. 

measured, — calculated 
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COO. 3 

' THETfl (DEG) 


Figure 4.15. Radar cross section pattern at f=ll . 95 GHz 

"i a 

E =0, <(>=0 plane. 

measured, - — calculated 




Figure 4.16. Radar cross section pattern at f=8.02 GHz, 

E 1 =<(> , <t>=90° plane. 

measured, — calculated 
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THETA (DEG) 


Figure 4.17. Radar cross section pattern at f=9.98 GHz 

E 1= <f>, <t>=90° plane. 

measured, — calculated 



RADAR CROSS SECTION (DBSM) 



THETA (DEG) 


Figure 4.18. Radar cross section pattern at f=ll . 95 GHz, 

E 1= <|>, <(»=90 o plane. 

measured, — calculated 
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Figure 4.19. Radar cross section pattern at f=8.02 GHz, 

E 1 =e, <(>=90° plane. 

measured, — calculated 
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RRDRR CROSS SECTION (DBSM) 



Figure 4.20. Radar cross section pattern at f=9.98 GHz, 

E 1 =0, 4>=90° plane. 

measured, — calculated 
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THETfl (DEG) 


Figure 4.21. Radar cross section pattern at f=11.95 GHz 

E 1 = 0 , 4>=90° plane. 

measured, — calculated 



Finally one notes that the rim scattering analysis performed here 
is valid for aspect angles away from the shadow boundaries of the 
diffracted field. For aspect angles close to the shadow boundaries of 
the edges (0 close to 90 degrees) the approach described here should be 
modified. This modification will not be given here. 

Frequency scans are done in the horizontal (<j>=0) plane, for the two 
different polarizations and the angle (9) fixed at four different 
values, namely, 9 equals 0, 15, 30 and 45 degrees. The measured (solid 
line) and calculated (dashed line) radar cross section results are 
plotted on the same graph and given in dBSM. The phase variation of the 
radar cross section is also shown in the figures. In order to be able 
to make a phase comparison, the measured data has been processed so that 
the two results have the same phase center. The frequency spectra are 
then inverse Fourier transformed to obtain the time domain responses. 

As explained before the time domain response is not an impulse response 
although it has been denoted as such on the plot for descriptive 
purposes Therefore, in order to make the comparison easier, the curves 
corresponding to measured and calculated time domain returns are shifted 
by an equal amount from the center line. The scale in the time domain 
is dimensionless and should be taken as a relative scale. Finally, due 
to the malfunctions of electronic instrument during measurements, there 
are some glitches in the measured data. They are kept as they appeared 
originally, since their presence does not affect the overall 
characteristics of the curves. 

As in the case of angular patterns, the amplitude and phase of the 
measured and calculated RCS results as a function of the frequency are 
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in good agreement especially for a 4> polarized field (Figures 4.22 
through 4.29). Also, it can be seen from the time domain plots that 
almost all of the measured return is from the rim at the open end. The 
absorber put at the back very effectively removed any internal cavity 
effects. 

* 

In the case of horizontal (e) polarization, (Figures 4.30 through 
4.37) the time domain plots show a relatively large return occuring later 
in time than the return from the open end. The variation of the position 
of this return with different values of (9) implies that it is coming 
from the external surfaces of the cavity model. Therefore as shown in 
Figures 4.32 through 4.37, the comparison between calculated and measured 
radar cross section results is not good. This also explains the 
discrepancy in the angular patterns of Figures 4.13 through 4.15. In 
order to compare the returns coming from the open end only, the unwanted 

return was gated out from the measured time domain data as indicated in 
each figure. The actual measured and calculated returns are shown with 
constant shifts, as before. These time domain plots were then converted 
into the frequency domain and shown on the same graph. As can be seen in 
Figures 4.32, 4.34, and 4.36 the calculated (long-dashed line) and 
gated-out measured (short -dashed line) results agree fairly well since 
they both correspond to the returns pertaining to the rim of the cavity. 

As explained before, the time domain responses shown are not the 
actual impulse responses because of the bandlimited nature of the data 
in frequency domain. If one had more frequency domain information then 
the time domain results would be closer to the true impulse responses. 
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radar cross SECTION (DBSM) 

-L|0. -30. -20. - 10 . o. 1C 


FREQUENCY IN GHZ 



FREQUENCY IN GHZ 


Figure 4.22. Variation of radar cross section with frequency 

E 1 =<j>, $=0, 0=0. 

measured, — calculated 
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Figure 4.25. Inverse Fourier transforms (i.e., time domain plots) of 
the results in Figure 4.24. 
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IMPULSE RESPONSE 



RADAR CROSS SECTION (DBSM) 


FREQUENCY IN GHZ 


1 



FREQUENCY IN GHZ 

Figure 4.28. Variation of radar cross section with frequency 

E 1 =J, 4>=0 » 9=45° 

measured, — calculated 



MPULSE RESPONSE 



Figure 4.29. Inverse Fourier transforms (i.e., time domain plots) of 
the results in Figure 4.28. 
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Figure 4.30. Variation of radar cross section with frequency 

E^e, <(>=0, 0=0° 

measured. — calculated 
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Figure 4.31. Inverse Fourier transforms (i.e., time domain plots) of 
the results in Figure 4.30. 






Figure 4.33. Inverse Fourier transforms (i.e., time domain plots) of 
the results in Figure 4.32. 
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Figure 4.35. 


Inverse Fourier transforms (i.e., time domain plots) of 
the results in Figure 4.34. 
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IMPULSE RESPONSE 
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Figure 
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Figure 4.39. Inverse Fourier transforms (i.e., time domain plots) of 
the results in Figure 4.38. 
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In order to show this, calculations corresponding to Figure 4.28 were 
repeated and expressions of Equation (3.15) were extended down to 10 
MHz. The resulting frequency domain plot is shown in Figure 4.38. 

Then, this data is inverse Fourier transformed, and the time domain 
result is shown in Figure 4.39. This result is closer to an actual 
impulse response and clearly shows the single, double and triple 
diffractions from the vertical edges of the open end. Note that the 
double order diffractions originating from both edges return to the 
receiver at the same time. 

ii) Results for the Interior Cavity Effects 

As described at the beginning, the coupling of incident energy into 
the interior waveguide modes and their subsequent radiation after 
undergoing multiple reflections between the back wall termination and 
the open end comprise the cavity effects. 

Since it is difficult to experimentally isolate the open end rim 
scattering from the cavity effects, the calculations include both 
effects for comparison purposes; namely, the results of (3.15), (3.134), 
(3.135), (3.140) and (3.141). The measurements and calculations were 
performed in two categories as in Chapter III; namely, aspect angle and 
frequency scans with the results presented in dBSM. 

The waveguide model is large in terms of the wavelength; therefore, 
a large number of propagating modes can exist inside the cavity. For 
example, at the frequency of 10 GHz, there are 152 propagating modes (86 
TE nm , 66 7W nm ) as well as an infinite number of evanescent modes. Since 
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the waveguide axial length is very long with respect to the wavelength, 
the evanescent modes will decay very significantly such that they can be 
neglected. Then the next question is whether the propagating modes have 
a preferred direction of radiation in that one can choose a few special 
ones rather than all the propagating modes and include only those in the 
analysis. It is has already been mentioned in the previous chapter that 
modal radiation is dominant along each of the discrete modal ray plane 
wave directions of that mode. Since each mode has discrete radiation 
directions, some will radiate strongly close to a desired direction, and 
others will not. This being the case, one can anticipate that only a 
few modes are significantly excited by the incident plane wave which 
radiate significantly in the backscatter direction. This claim has been 
checked numerically, and the results are shown in Figure 4.40. An 
open-ended rectangular waveguide with dimensions equal to the 
experimental model used in this study is analyzed at 10 GHz. The 
backscatter field is calculated in (x-z) plane by varying the aspect 
angle (0). Note that only the modal effects are included and the 

a 

incident field is assumed to be ^-polarized. It is found that only the 

TE modes are excited in this plane for this polarization. Their plane 
no 

wave directions or modal ray angles are tabulated in the figure. The 
scattered field is calculated by including all modal contributions as 
indicated by the solid line. For comparison purposes the contributions 
of TE no (n=l,...5) modes and TE nQ (n=6...13) modes are shown in the same 
figure. Note that the first five modes radiate strongly and almost 
replicate the solid curve in the region close to their modal ray angle 
directions. The remaining seven modes are major contributors for large 
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MODAL RAY ANGLES OF TE no MODES 
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Figure 4.40. Importance of modes whose modal ray angles are near the 
angle of incidence. 

•l A 

f=10 GHz, E =<{>, 4>=0 plane. 





aspect angles. After observing this fact, the research focused on 
finding a rule of thumb procedure which could be used to select the 
minimum number of modes needed for a given direction of incident plane 
wave. From numerical results and comparisons with the experimental 
data, it was found that in the horizontal ( 4>=0) and vertical ( <|>=90 o ) 
planes only three modes were sufficient for the cavity under test. 

These three modes are selected based on their modal radiation direction 
such that they are closest to the incident plane wave direction. The 
aspect angle scan results are shown in Figures 4.41 through 4.64. As 
done in the previous section, the experimental results are indicated by 
the solid line and calculations by the dashed one. In each case, a 
comparison is provided for the calculations which include all 
propagating modes versus three propagating modes. In all cases, a 1 
dB/bounce energy loss is assumed to model the imperfection of the 
conductivity of the model used in the measurements. 

Since this new concept proved to be so valuable in the principal 

planes; namely, <j>=0 and <(>=90° planes, it is next applied to <j>=45° plane 

to see if it fails when the incident plane wave direction is not aligned 
with the structural symmetry. The aspect angle scan is calculated in 

A A 

the <j>=45° plane for both the <(» and 9 polarized incident fields. The 
results are shown in Figures 4.65 and 4.66. In this case, the 18 
preselected modes are compared with the complete 152 modes as shown in 
each figure. Note that more terms are necessary in this general case as 
indicated by the results shown in Figure 4.67 where the six mode result 

is compared with the 152 mode one. Even so it is clear that one can use 

far fewer modes than the complete propagating mode set. 
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Figure 4. 



. Radar cross section pattern at f=8.00 GHz, 
A 

E =(j>, $=0 plane, 

measured, — calculated 
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Figure 4.42. Calculated RCS pattern corresponding to Figure 4.41 

All modes are included 

xxx Only 3 modes are included 





Figure 



.43. Radar cross section pattern at f=10.00 GHz, 
"1 * 

E =<|>, ((>=0 plane. 

measured, calculated 
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Figure 4.44. Calculated RCS pattern corresponding to Figure 4.43. 

All modes are included 

xxx Only 3 modes are included 
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Figure 4.45. Radar cross section pattern at f=12.00 GHz, 

a 

E =<f>, <J>=0 plane. 

measured, — calculated 
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Figure 4.46. Calculated RCS pattern corresponding to Figure 4.45 

All modes are included 

xxx Only 3 inodes are included 
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Figure 4.48. Calculated RCS pattern corresponding to Figure 4.47 

All modes are included 

xxx Only 3 modes are included 







Figure 4.50. Calculated RCS pattern corresponding to Figure 4.49. 

All modes are included 

xxx Only 3 modes are included 
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Figure 4.52. Calculated RCS pattern corresponding to Figure 4.51 

All modes are included 

xxx Only 3 modes are included 




Figure 



.53. Radar cross section pattern at f=8.02 GHz, 

E 1= 0, <}>=90 o plane. 

measured, — calculated 
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Figure 4.54. Calculated RCS pattern corresponding to Figure 4.53. 

All inodes are included 

xxx Only 3 modes are included 
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Figure 



.61. Radar cross section pattern at f=10.01 GHz, 

E 1= k <t>=90° plane. 

measured, — calculated 
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Figure 4.63. Radar cross section pattern at f=ll . 96 GHz, 

E 1 =<j>, «|»*90 o plane. 

measured, — calculated 
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Figure 4.64. Calculated RCS pattern corresponding to Figure 4.63. 

All inodes are included 

xxx Only 3 modes are included 
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Figure 4.65. Radar cross section pattern at f=10.00 GHz, 

-4 * 

E =<{>, <j>=45° plane. 

All 152 modes are included (86 TE, 66 TM modes) 

xxx Only 18 modes are included (9 TE, 9 TM modes) 






Figure 4.67. Radar cross section pattern at f=10.00 GHz, 

E 1= e, 4>=45° plane. 

All 152 modes are included (86 TE, 66 TM modes) 

xxx Only 6 modes are included (3 TE, 3 TM modes) 
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msforms (i.e., time domain plots) of 
ire 4.68. 
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Figure 4.71. Inverse Fourier transforms (i.e., time domain plots) of 
the results in Figure 4.70. 
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Figure 4.73. Inverse Fourier transforms (i.e., time domain plots) of 
the results in Figure 4.72. 
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Figure 4.74. Variation of RCS with frequency E 1 ^, <(>=0, 9=45°. 
measured —calculated 
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Figure 4.75. Inverse Fourier transforms (i.e., time domain plots) of 
the results in Figure 4.74. 
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Figure 4.79. Inverse Fourier transforms (i.e., time domain plots) of 
the results in Figure 4.78. 
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Figure 4.81. Inverse Fourier transforms (i.e., time domain plots) of 
the results in Figure 4.80. 
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Figure 4.83. Inverse Fourier transforms (i.e., time domain plots) of 
the results in Figure 4.82. 
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Figure 4.84. The separate returns of modes for the case 
in Figure 4.75. 



Finally, the frequency scan results are calculated and compared 
with measurements in Figures 4.68 through 4.83. As in the previous 
section, the time domain results are deliberately shifted to better 
illustrate the comparisons. Note that in each case the comparisons are 
very good. 

In Figure 4.84, the contribution of several modes is plotted 
separately for the geometry associated with the results of Figure 4.75. 
It is seen that as the mode number increases the given mode undergoes 
more bounces inside the waveguide, travels a longer distance, and 
therefore its return is received later in time; i.e., a stronger modal 
dispersion. Also note that the modal terms add up to form ripples in 
the total result from the cavity effect. If not understood correctly, 
one might think this is indicating different scattering centers, which 
would not be true. This shows that modal propagation should carefully 
be traced if one wishes to obtain the true response for the cavity 

A/»f r 

C I I CUb J • 

e) Electromagnetic backscattering from a circular waveguide 
cavity: 

The EM backscatter results based on (3.5), (3.6), (3.151) and 
(3.152) are calculated here. The rim and cavity effects are separately 
shown together with their superposition in Figures 4.86 through 4.89. 

The main contribution to the backscatter return is the cavity effect due 
to the perfectly-conducting termination; whereas, the rim scattering 
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Figure 4.85. Circular waveguide cavity terminated by a short circuit. 

contribution results in a perturbation to the cavity effects. The 
radius of the waveguide is taken as 3.34 wavelengths which allows 115 
modes to propagate. The Figures 4.86 and 4.88 illustrate the 
contributions of all propagating modes; whereas in figures 4.87 and 4.89 
only the modes are included with modal ray angles inside a 10° 
neighborhood of the observation direction. It is clear from these plots 
that by including only the few significant modes, one can substantially 
reduce the amount of calculations without seriously sacrificing 
accuracy. It is also noted that the length of the waveguide cavity is 
10 wavelengths; therefore, all evanescent mode contributions are 
neglected. 
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Figure 4.86. Backscattered field from a circular waveguide cavity, as a 

function of the angle from the axis. All propagating 

i A 

modes are included. E ff <f> . 
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Figure 4.87. Backscattered field from a circular waveguide cavity, as 
a function of the angle from the axis. Only modes with 
modal ray angles closer than 10 degrees are included. 
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Figure 4.89. Backscattered field from a circular waveguide cavity, as 
a function of the angle from the axis. Only modes with 
modal ray angles closer than 10 degrees are included. 



CHAPTER V 


CONCLUSIONS 

A relatively simple, efficient and sufficiently accurate analysis 
of electromagnetic scattering by a class of electrically large, 
open-ended waveguide cavities was developed in this work using a hybrid 
combination of high frequency, modal and multiple scattering methods. 

The cavities are composed of waveguide sections in which the Helmholtz's 
equation is separable; therefore, the EM field inside each section can 
be written as a sum of waveguide modes. These modal waveguide field 
expressions are expressed in terms of "modal rays" through asymptotic 
approximations to the modal functions. The scattering properties of the 
discontinuities formed by the junctions between the sections are 
analyzed using high frequency techniques together with the modal rays. 
The latter techniques employ the GTD and the ECM which require a 
knowledge of the pertinent diffraction coefficients that are available 
from the asymptotic solutions to appropriate canonical problems. In 
some situations the analysis also employs high frequency approximations 
based on the PTD and its modifications. These high frequency techniques 
used in conjunction with the "modal rays" lead to a relatively simple 
description of the isolated junction scattering matrices. These 
individual junction scattering matrices are then combined in a 
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self-consistent multiple scattering approach to get the total scattering 
effect of the cavity. 

Basically, the scattering matrices are analyzed in two categories. 
The reflection type scattering matrices include an integration of the 
GTD based equivalent currents over the edges of the junction apertures. 
It is important to note that the previous ray-optical techniques 
[10,14,20] to treat the reflection from the open-end of waveguides 
employed the GTD based equivalent electric and magnetic line sources in 
free-space. That procedure is applicable only to those geometries where 
one can sum up the radiation from equivalent sources and their images, 
such as parallel -plate and circular waveguide geometries. In the 
present work, the equivalent sources radiate in the presence of 
waveguide walls and therefore they are magnetic line sources and 
magnetic line dipoles. Hence, this proposed procedure can be 
generalized to many other goemetries which lack the symmetry properties 
required by the previous methods. The transmission type equivalent 
currents are computed via a Ki rchhoff-Huygens 1 approximation to the 
aperture field. It is shown in Appendix H that this aperture 
integration can be corrected in the same way as the PO is corrected via 
the PTD. However, the integration process does not give any physical 
insight to the problem, and as the number of propagating modes increases 
with frequency, it becomes cumbersome and inefficient. 

In this research, new approaches to substantially improve the 
efficiency of the above mentioned aperture integrals were investigated. 
It was determined that for the open-ended rectangular and circular 
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waveguides the modal radiation from the open end is strongest along 
their modal ray directions. Therefore, for a given radiation direction 
one can include only the radiation of those modes whose modal ray 
directions are closest to the radiation direction. In Section 3.2, it 
was shown that for the cavities formed by parallel -plate waveguides, if 
one includes only the contributions of modes with modal ray angles which 
are less than a fixed angular distance from a given observation 
direction, then the error in the backscattered field stays bounded with 
increasing frequency, even though the number of propagating modes 
increases with frequency. In Chapter IV, numerical and experimental 
results were compared to assess the accuracy of the analysis, and to 
indicate that the contribution of only a few modes with modal ray 
directions closest to the observation direction accurately approximates 
the contribution of all the propagating modes. This is an important 
result, because it combines modal radiation with ray-optics, and 

tharofnro it ran ho annli pH to manv Hiffprpnt anrt rnmnlpy wavpnui rlf* 

»..»■ wi « * w ~ >- r r * • — — v — r — - - ■ - — - - ^ - 

geometries to effectively select only the few significant modes from the 
entire set of propagating modes. 

The accuracy of the scattering from waveguide cavities is 
determined in turn by the accuracy of the scattering matrices involved. 
In this research, the accuracy of the analysis for scattering matrices 
has been verified by comparison with other analytical and experimental 
results on certain cavity geometries. It is noted that there were very 
few dependable numerical and experimental results available for the 3-D 
cavities formed by sections of linearly tapered and uniformly curved 
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waveguides even though they are very often encountered in practice. The 
experimental model with linearly tapered walls which was treated in 
Chapter IV had a small interior discontinuity, and, as expected, it did 
not seem to significantly influence the final results. As a part of 
future research, practical cavity models with more pronounced interior 
discontinuities will be built and analyzed. 
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APPENDIX A 


MODAL FIELD EXPRESSIONS IN A RECTANGULAR WAVEGUIDE 



Figure A.l. Rectangular waveguide geometry. 


As described in [44], the modes in a rectangular waveguide as shown 
in Figure A.l can be classified into sets of fields transverse to a 
coordinate direction. 

A 

The mode sets transverse to the axial (z) direction are important 
and widely used, because it applies to uniform nonrectangular 
cross-section guides. However, in many problems, mode sets transverse 

A A 

to x or y coordinate may be more suitable. 
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These sets can be calculated from an axial or transverse vector 
potential [44]. The TE nm and TM nm to z type modal field expressions are 
given by: 

TEnm TM nm 
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n a +m b 


h = N Y 
z nm nm 


7“ cos n x cos m. y 
J^nm a b 
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(A.la;b) 
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Note that for both modes: 
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where the normalization factor N is defined such that 

nm 

a b 

1 ★ A 

7 / / 5 x fi • z dy dx = 1 , 

L Q J ' nm nm * 


(A. 13) 


where e and h are the tangential components of the field as follows, 

nm nm * 

+j B n z A +j Bz 

E . . = [e x+e y + e z"| e = [e ±e , z] e nm (A. 14) 
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and denotes complex conjugation. 

a 

The upper or lower sign represents a mode propagating in (+z) or 

A 

(-z) direction, respectively. 

The modal field expressions in the rectangular waveguide can be 

A 

decomposed into four plane waves. The transverse (to z) components of 

A 

the fields propagating in +z direction are written as follows: 
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where 


u = 


m b for TE modes 
n a for TM modes 


(A. 18) 


and 


v = 


n a for TE modes 
-mb for TM modes 


(A. 19) 


The modal sets transverse to x-axis can be written as follows: 
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From (A.l) through (A. 25), one can easily conclude that 
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APPENDIX B 


MODAL FIELD EXPRESSIONS IN A PARALLEL PLATE WAVEGUIDE 

The modal field sets in a parallel plate waveguide of width a, shown 
in Figure B.l, can be obtained from (A.l) through (A. 12), by discarding 

A 

the field variation in y-direction. They are expressed as follows: 


x 



0 


Figure B.l. Parallel -pi ate waveguide geometry. 
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The normalization fabtor N now becomes 
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The rest of the parameters are obtained from (A. 9) through (A. 12), 
by letting m=0. 

Each of these modes can be decomposed into two ray optical fields 
(plane waves) by writing 
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(B.9) 


in the TE^ case, and 


H = y 


Nn 

r 




+e 


■jn a x 


■j v 


in the TM^ case. 

The propagation directions of these plane waves are called the modal 
ray directions and the angles between the modal ray directions and 
waveguide walls are the modal ray angles. 



Figure B.2. Ray picture of n^ mode. 
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APPENDIX C 


MODAL FIELD EXPRESSIONS IN A SECTORAL WAVEGUIDE 


y 



Figure C.l. Geometry of a sectoral waveguide. 


In this problem, it is possible to find mode sets transverse to the 
x-axis. The fields can be computed from the scalar function 


32 32 

nw 


3y 2 + 3Z 2 + 

* 2 - C ^ 2 






(C.l) 


with the appropriate boundary conditions. is the x-component of the 
magnetic vector potential A in the TM X case and the x-component of the 
electric vector potential F in the TE X case [44], 
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Case 1 


TE modes: These modes satisfy the relation E =0 and can be 

X x 

— A 

derived from F=xi|>. 


field components 

are given by; 
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The appropriate scalar function <|> is given by 
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The explicit expressions for the outgoing field components are given by; 
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E x * 0 (C.9) 
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For small tapering angle and large p, the Debye Asymptotic 
approximation for the Hankel function [45] 
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can be used. For small 7 ; cos - 1 7 + 7 , so the above relation becomes 
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which goes to the large argument form [44] for z»v. 
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If one makes use of (C.14) and (C.15) in the field expressions of 
(C.7) through (C.12), one obtains the following expressions (up to a 
common factor) for an outgoing wave: 
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By comparing (C.16) through (C.22) with (B.20) through (B .25) , one can 
see that the expressions of slightly tapered waveguide fields may be 
approximated by TE X , TM x -type rectangular waveguide mode sets. To do 
this, one has to make the following approximations: 
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together with (C.25) where the p-dependence cannot be approximated by 
(C.28), since it is the function determining the phase advance which is 
more sensitive to the approximations. 

Case 2 TM modes: These modes satisfy the relation H =0 and can be 
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derived from A=xtp. 
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The field components are given by: 
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The explicit expressions for the outgoing field components are given 


by: 
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Again using Debye's Asymptotic form for the Hankel function, one 
the following field expressions: 
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These modal expressions under the approximations defined by 
(C.26) through (C.28) become equivalent to the TM X modes in a 
rectangular waveguide. 

The ray picture of the modes can be obtained using a similar 
procedure performed for a rectangular waveguide. Each outgoing (or 
ingoing) mode corresponds to four rays which follow zig-zag paths inside 
the waveguide. 

In the TE X case, the ray optical fields are obtained from 
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In the TM X case, the ray optial fields are obtained from 


H p + H . J = 
p ♦ T 


L n r 
nm o 

4/ e n mP 


nw mir 

JCr*) (*;♦) -i s nm p 


mir 

_p<|> 0 + e nm 


+e 


mr mir 

jto) -j to) -j*p »P 


nnr 




-Jto) + J 

+e 


(♦»♦) ~^nm p 


mir 

P«l> 0 + B nm 


,mr 


+e 


jto) -j 


/ mir A 

(♦O*) ~j P l 


nm 


mir 

p 4> 0 - 0 nm 


(C.45) 


249 



If one looks at the projection of outgoing rays corresponding to 
n** 1 mode in the y-z plane, one obtains rays following zig zag paths 
along the guide bouncing from the tapered walls (via the approximation 
in (C.13)). The ray trajectories are tangent to the circular modal ray 
caustic whose radius is determined by the mode index as shown in Figure 
C.2. It is also noted that inside the circular modal ray caustic, the 
mode is cut-off and the ray representation is not valid. 

The projection of the ray picture onto the x-z plane is shown in 
Figure C.3, where the rays are bouncing from the parallel walls. 



sine 
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kp<l>o 


nm 
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( k ) 2 


Figure C.2. Projection of ray picture into y-z plane. 



Figure C.3. Projection of ray picture into x-z plane. 
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APPENDIX D 


MODAL FIELD EXPRESSIONS IN A 2-D LINEARLY TAPERED WAVEGUIDE 

The geometry pertaining to the linearly tapered waveguide is shown 
in Figure D.l. The expressions for the modal fields can be deduced from 
the sectoral waveguide by taking the width "a" to infinity; they are 
given as follows: 
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0 

(D.2a;b) 

L 

— — . -J^mP 

^ S1n v" e 

(D.3a;b) 

mir — [E_ jY o mir -J0mP 

p+o cos <f. 0 <t> e 

(D.4a;b) 
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(D.5a;b) 
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Figure D.l. Geometry of a linearly tapered waveguide. 


where 


L 

m 



‘t'o 


is the normalization constant with 


e om 


= { 


2 ; 

1 ; 


n=0 

n*0 


(D.7) 


(D.8) 


and 


m 



(D.9) 
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The next step is to write the modal field in ray-optical form. In 
the TE X case one can write (D.6a) as follows: 


E 

O 

>- 

j 

m 

. d> - Bp 

-j 

irm 

* = 2/0^ 

e 

9o Y m 

+e 

?o m 


Equation (D.10) shows that a modal field is composed of two ray 
optical fields. The ray picture of the modal fields is illustrated 
below [18], 



Figure D.2. Ray picture of n^ mode. 
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The ray field is obtained in the TMy case by writ 
f ol 1 ows : 


E 


x 


Y I — 

om . , x . nwr x 

=T J " e n/0 " J Uo* + e m p ) 

2j/ 3 m p _ e -e 
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APPENDIX E 


CIRCUMFERENTIALLY PROPAGATING MODES IN AN ANNULAR REGION 


y 



Figure E.l. Geometry of the annular region. 


In this appendix, circumferentially propagating modes in an annular 
region will be obtained using Green's function techniques. 

The direct determination of the Green's function for 
circumferentially propagating waves can be done in a way similar to that 
presented by Wasylkiwskyj [46] for the interior problem of a single 
circular shell. The Green's function is the solution of 


A _ a _a 1 a 2 , 2 . 

Ip 3p P 3p + p 7 34 z " + * > 


G(p,p' ) 


Mp-p'HU-V) 

“ ’ — | — 


(E.l) 
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with the boundary conditions 

G = 0 at p = a,b . (E.2) 

for the TM Z case. The TE Z case analysis can be performed similarly. In 
addition to (E.2), the radiation condition at has to be satisfied. 
The solution of (E.l) can be written as 

G(p.p’) = I * n (p)* n (p') (E.3) 


where g < j ) (4 , *4>* i v n ) satisfies the differential equation: 


d 2 2 " 

W + v 


g^ = -«(♦-♦') 


(E.4) 


with 


3 <l> 


as <|> -*• ± 


(E.5) 


and * (p) are the normalized eigenfunctions in the radial direction 
which satisfy 


p dp (Pdp) + k2 P 2 


i|»(p) = v 2 \|/( p) 


(E.6) 


with 


i|> = 0 at p=a,b 


(E.7) 


The solution of (E4) and (E.5) is given by [46] 
9* (+»*'; v ) = — 2 JF 0 


(E.8) 
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and ^(p) can be found from the radial Green's function gp(p,p';v 2 ) 


"1 d d 

p 3p (P3p) + k2 ) - y 

P 2 " 


9 p = ~ s (p-p' ) 

P 


(E.9) 


with 


g = 0 at p=a,b and y=v 2 
p 


(E. 10) 


<|> n (p) are obtained from g p via the completeness relationship 


p'«(p-p') - E * n (p)* n (p') = / g (p,p';y)dy 


(E.ll) 


n 


where C is the contour encircling the singularities of g . 

p 


The solution to (E.9) and (E.10) is given by: 


g p (p,p';v 2 ) = +j } 

[H^ 1) (kp < )H^ 2) (kb)-H^ 1) (kb)H^ 2) (kp < )][H^ 1) (kp > )H^ 2) ( ka ) -H^ 1 } ( ka )H^ 2 } ( k P> ) ] 

[H^rkalH^fkhl-H^ fkblH^ (kail 

L V ' V V' 'v' ,J 

(E.12) 

The expression in (E12) can be written in a more suitable from by 
letting 

[Hi 1) (Kp^)H^ 2) (kb)-H^ 1) (kb)Hj 2) (kp )] sf (p ) (E.13) 


so that 


9 p ( p»p‘ ; v 2 ) 


j. 

4 H^(kb) 

V 


f(p<)f(p>) 

f(a) 


(E.14) 
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Using this expression in (Ell), 




J *n(P)*n( p,) = + l 2 V n h (2) 


H vn )(k3) f(p)f(p') 


H. ; (kb) 


WT 

3v I v=v n 


(E.15) 


Note that dp=2vdv is used in obtaining the residue sum in (E.15), where 
residues are evaluated at the isolated, simple poles [47] given by 

f (a) = 0 . (E-16) 

The Green's function of (E.3) then is written as 


(2) 

, -jv U-*'| \ < ka ’ f(p)f(p‘) 

G(p.p') = T. j e .,(2),.. . »f(a) 


H* ' (kb) 


(E.17) 


3v I v= v n 


This result can also be obtained using an application of Watson's 
transformation [48,49] on 2ir-periodic eigenfunctions. This will also be 
illustrated for the sake of completeness. 

The Green's function of (E.3), can be written as follows: 


G(p,p') = Z • ( ♦)* ( *')g (p,p';X m ) 

m p 


(E.18) 


where eigenfunctions in the 4> direction satisfy the boundary condition 
of being periodic in 0<<f><2ir. These eigenfunctions are given by 


fc m 

*„,(♦) = W cos '"('l'-*') 


where 




1; m = 0 
2; m * 0 


(E. 19) 


(E.20) 


9 p is given by Equation (E.12). 
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Substituting (E.14) and (E.19) in (E.18) 

H ' 2) (ka) 




m=-co 


m (kb) 


f(p)f(p') 

— TUT~ 


v=m 


(E.21) 


This summation can be written as the following integral 


G(p,p') = ‘2T / I(v) e 

C=C + +C_ 


jv(<M>') e Jwr 


sinvir 


dv 


(E.22) 


where C encloses zeros of sinvir in the counterclockwise direction. In 
(E.22) I(v) is given by: 


.(2) 


I(v) " “8J h (2) 


1 H v (ka) f(p)f(p') 


H^'(kb) 


w 


(E.23) 


and has no singularities in C, as shown in Figure E.2. 


Im v 



Figure E.2. Contours used in the Watson transformation. 
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In the figure, the singularities of I(v) are assumed to be encircled by 
the contour C' (C'=C' + +C' ) in the complex v-plane. Deforming the 
original contour C into C 1 and evaluating the integral as a residue sum, 
one obtains 


tt 00 -jv n U-<|>' ) 

G( p,p' ) « E e 
n=l 


j V 

e 

sin v w 


H^ 2) (ka) 
v n 

H (2) (kb) 

V 

n 


f(p)f(p‘) 

9f(a) 

3v v=v n 


tt 00 ~j v n ( 1 ) 

+ 8J S e 
n=l 



TT 


sin v ^ 
-n 


Tt 


H^ 2) (ka) 
v -n 

(kb) 


f(p)f(p') 

9f(a) | 

3v I v=v n 

(E.24) 


The first sum corresponds to upper half plane poles which represent 

A 

modes propagating in +()» direction and the second summation is obtained 
from the residues of the poles in the lower half plane representing 

A 

modes propagating in -<|> direction. By utilizing the symmetry properties 
of Hankel functions both kinds of modes can be included in one summation 
as follows: 


G(p,p' ) 


TT 

8J 


l 

n=l 


-j ^ I 

e 



sin Vp tt 


f(p)f(p') 
3f(a) I 
3v I v=v n 


(E.25) 


Also using 


Jv n 11 

e °» -2jv irp 

= 2j E e n 


Sin Vpir 


(E.26) 


p=0 


one notices that each term in (E.25) corresponds to an infinite number 
of waves reaching the observation point after p complete encirclements. 
Therefore, this representation of Green's function contains multiply 
encircling waves, which are included in the 2ir-periodic eigenfunctions. 
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The first encirclement (p*0) gives the following expression for the 


Green's functions; 

( 2 ) 

, - -JvJt-VI \ (ka) f(p)ftp') 

G(p»p‘) - T S . e (?\ 3f (a) | .27 ) 

n= i (kb) 3 V | v=v n 


which agrees with (E.17) above. 

In deforming the contour C into C' , one has to show that the 
integrand of (E22) is exponentially small on the semicircular contours 
E+ shown in Figure E.2 at infinity. This can be done using the Debye 
Approximations to Hankel functions [45,49] and the fact that 


lim 


e jv I*-.). | 

sin v<j> 


M e -R(*|sine| + (ir-e)sin9) 


(E.28) 


i d 

where M is bounded and v=Re . 

As is observed in (E.17) and (E.27), the eigenvalues of 
circumferentially propagating modes are determined through the zeros of 
f(d), v/hich is given by 

f (a) = H (1) (ka) H (2) (kb) - H (1) (kb) H (2) (ka) . (E.29) 

v n v n v n v n 


Using asymptotic approximations to Hankel functions, the zeros of 
f(a) can be located in the complex v-plane. As shown in Figure E.3, the 
zeros are symmetrically located on the real and imaginary axes for real 
values of the wavenumber k. The purely imaginary eigenvalues correspond 
to evanescent modes which decay with the distance in the circumferential 
direction. The real eigenvalues correspond to propagating modes, which 
travel inside the annular region. The magnitudes of these real 
eigenvalues are less than the electrical radius of the outer shell (ka). 
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Figure E.3. Location of eigenvalues in the complex v-plane. 


From (E.17) and (E.27), the expressions for the eigenfunctions 


corresponding to the eigenvalues are given by; 


Q n = A n (ka.kb) e 




Where A® is a normalization factor determined by the inner and outer 
n 

radii. To get the ray picture of propagating modes, one needs to look 
at the eigenfunctions corresponding to real eigenvalues. If |v n |<kb<ka, 
then the Hankel functions in f(p) of (E.30) can be approximated by their 
Debye Asymptotic forms given in (C.13). This approximation gives 
incoming and outgoing rays, bouncing from the inner and outer shells and 
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staying tangent to a circle with electrical radius v p as shown in Figure 
E.4. The modal ray angles are given by the following expressing: 


6 


c 

n 


TT_ 

2 


-sin 



upper shell 
_lower shell 


(E.31) 


This ray picture is similar to the ray picture in a parallel plate 
waveguide in the sense that the rays bounce from both walls. Therefore, 
these modes will be referred to as regular modes (RM). 

For kb<|vjka, one can use the following spatial filtering property 
of Bessel functions; 

J y (x) *» 0 for |v|>x (E.32) 


MODAL RAY BOUNCING 
FROM BOTH WALLS 



Figure E.4. Ray picture of regular modes. 
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in the expression for eigenfunctions. Then (E.30) can be written as; 

-J' v n UI 

Q n 3 A n ( ka » kb ) e (-J2) N V ( kb )J v ( k p) (E.33) 

n n 

where N Vn (x) is the Neumann function of order v n « By use of (E.32), the 
mode expression in (E.33) reveals that most of the energy in this case 
is attached to a region close to the outer shell. The ray picture of 
these modes is shown in Figures E.5 and the caustic circle is in the 
annular region, resulting in the ray bounce occuring only on the outer 
shell. This ray interpretation is also discussed in [50] and similar to 
the ray interpretation of whispering gallery modes inside a circular 
cylindrical region [46,51]. Therefore, these modes will be referred as 
whispering gallery (WG) modes. The modal ray angles of these WG modes 
are given as follows: 

w * _i , v n, 

6 n = 2 -sin 1 (i^-) . (E • 34) 

MODAL RAY FOR A 
WHISPERING GALLERY MODE 



RAY CAUSTIC 

Figure E.5. Ray picture of whispering gallery modes. 
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APPENDIX F 


MODAL FIELD EXPRESSIONS IN A CIRCULAR WAVEGUIDE 


The circular waveguide has a radius a, as shown in Figure F.l. 
The field sets can be written as TE Z and TM Z types. The propagation 
constants of the TE Z and TM Z modes are given by 


0 nm 


nm 2 

k 2 - (— ) 


and 



K nm 2 

k 2 - (— ) 


(F.l) 


(F.2) 


where p is the m tn root of the n tn order Bessel function, namely 
nm 


J„(p ) = 0 
n VK ntrr 


(F.3) 


* t h 

and p is the m root of the derivative of the Bessel function of 
r nm 

order n; as follows: 


j' (p' ) = 0 
n VK nrr 


(F.4) 


265 



X 



Figure F.l. The circular waveguide geometry. 


The field expressions are given as follows apart from the phase factor 

A 

along the z direction: 
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TE modes 


TM modes 


h 


z 




rcos n<|> 
lsin n<j> 


0 


(F.5a ;b) 


e 

z 


0 


jN nm (Pn!H) 2 J n (fnrnP) 
nnr a ’ n ^ a J 


rCOS n<j> 
lsin n<j> 


(F.6a;b) 


i i , 


h 

p 

N 1 ^nmPnm l c ? s n< ^ 

nm a nl a J lsin n * 

. e 4> 
^e,nm 

(F.7a;b) 

% 

1 1 

n 1 n ^nm j fPnmP'i r-sin n<j> 
nm p V a J l cos n+ 

e P 

^e.nrn 

(F. 8 a;b) 

e 

p 

^h,nm% 

N 0 nm p nm ,vf p nmP'| f cos n + 
nm a V a > l sin n<(> 

(F.9a;b) 

% 

-Z. h 

h,nm p 

N n ®nm ,i rPnmP'i r-sin n<|> 
nm p V a J 1 cos n * 

(F.10a;b) 

1 ■ 

'"h.nm 

If -7 
nm 


(F.llajb) 

Z 

e,nm 


®nm 7 

nr ° 

(F.12a;b) 


2 


(F.13) 

J n( p nnJ V^“ w ^nm e on( p nm^ 

ITT 

N 

nm 


2 

. (F.14) 


P nm ^n ( p nnJ V^^nmSn 
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The normalization factors N and N yield unit power in each mode. 

nm nm 

Here, e Qn = 2 for n = 0, and 1 for n > 0. 

To obtain a ray picture for the typical modal field whose p and z 
variation is given by 


J ( E, p) e 
n v nnr' 


-j Y | Z | 

nm 1 1 


(F. 15) 


where 


a«5 


nm 


and 


p for TE modes 

nm z 

p for TM modes 

nm z 


(F.16) 


nm 


8 for TM modes 

nm z 

6 for TM modes 

nm z 


(F. 17) 


one may decompose J p (x) as 


«' 2, (x) + h! 1) (x) 


J n<*> * 


(F. 18) 


Thus, the expression in (F.15) becomes 

\ H I 2) (5 nm p) e JYnm|Zl + \ H^ 1 ^ ( 5„ m p) e JYnm ' Z ' 


nm r 


’nm r 


(F.19) 


From the large argument asymptotic form of the Hankel functions, it can 
be determined that the first term represents conical rays propagating 
away from the axis of the waveguide and the second term represents 
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conical rays propagating toward the axis of the cylinder. The described 
ray picture has caustic along the guide center, therefore it is valid 
away form it. The modal ray picture is sketched in Figure F.2 below. 

I 

The modal ray angles . « nm and 6 nm are as follows: 
i i p nm 

(— ) (F - 19 > 

and 

i p nm 

s„ m - sin' 1 (— ) . (F.20) 

nm v 

It is also noted that the modal field is not completely transverse 
to the conical ray propagation direction, since the radial component of 
the field has contribution along the conical ray propagation direction. 
However, for the modes sufficiently away form cut-off this contribution 
is small and can be neglected to obtain a ray optical representation. 


DIVERGING RAY CONES 



Figure F.2. Conversion of waveguide modal field into the conical ray 
field. 
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APPENDIX G 


WAVEGUIDE EXCITATION PROBLEM 


r 



Figure G.l. Waveguide Geometry 

Let a source S be located inside a waveguide as shown in Figure 
G.l. S can be an electric current source J or a magnetic current source 
M, which generates E 4 , H 4 . Thus, 

r, t = ea; (i 4 h 4 ) (g.i) 

p r K K 

*■ A 

are the fields generated by S in the ±n directions. In this summation p 
is the compact summation index representing the double summation over 
mode indices V and "m". It is of interest to find (E 4 ,H 4 ). Here, E 4 
will be determined via an application of the reciprocity theorem as in 
[ 1 ]. 
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The modal fields can be decomposed into transverse and axial 
components as follows: 


-+ 


H~ = (±H . + H ) e p 
p v pt pn 


(G.2a) 


-+ +JM 

E- = (E . ± E ) e p 
p pt pn 


(G.2b) 


Case (1): It is of interest to find the strength of the mode E”, 
hf generated by S. (Here q is the compact index representing the mode 
indices "ij"). 

From the reciprocity theorem: 


( JJJ E + «J dv for S=J 
1 V q 

-HI H + «M dv for S=M 
[ V q 


II [ExH -E xH]*n ds = 
S.+S +S q q 


+ - w 


(G.3) 


where E. H = E + , H + on S. and E, H = E“, H“ on S . Also, E + x n ] = 0. 

M s' 
w 

Thus, the above reciprocity relation simplifies to: 


// [i + xHq-E*xH + ] - n ds + ;/ [rxii^xH-] • (-S) ds 

s s 


III E + * J dv for S=0 

V q . (G.4) 

-HI H + *M dv for S=M 

V q 
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Substituting E , H and E , H~from (G.l) and using the orthogonality 


condition: 


-+ -+ 


|J [E- <»;]•« ds 


= 0 for q * p 


one obtains a relationship in terms of only the qt* 1 mode (i.e., p=q 
case) 

, . . _ - - - -j 2e n n + 

If A [ (E +E )x(H .xH )e q -(E +E )x(H xH )e q ] 
q L qt qn qt qn qt qn qt qn J 

+ 

-// A“[(E .-E )x(H .+H )-(E .+E )x(-H .+H„ )]«n ds 

q 1 '' qt qn v qt qn qt qn qt qn J 


or 


ffl E O dv for S=J 

V q 

-J J J H + .M dv for S=M 

V q 


-2 A" If E . xH »n ds = 
n qt qt 


III e J8q " dv for 5=J 

. _ _ -J M 

-Iff ( H 4 - +h )*M e q dv for S=M 

\l qt qn 


Thus 


A" = - 


q 2 If E xH -n ds 
s qt qt 


- _ _ -j 3 n 

III (E .+E )*J e q dv for S=J 

V & ^ 

- - . -je n n 

■fll (H +H )*M e q dv for S=M 
\ j qt qn 


where A” is the excitation coefficient of the mode E" traveling to 

q q 

right of the source, i.e., in the region n<n shown in Figure G.l. 


(G.5) 


•n ds 


(G.6) 


(G.7) 


(G.8) 


the 


Case (ii): If E~, H“ is chosen as the waveguide mode of interest 

q q _ . 

as in Figure G.2, then the coefficient A* of. the mode E* traveling to 

q q 

the left of the source (n>n + ) is given by: 


1 




n ds 


- - -J Bn - 

^W Je dv 

-iii <-v+v- s e ' J8q, ' dv 


(6.9) 


Note that fj E qt xH qt .nds=//E qt xH qt .nds=//E qt xH qt .nds where S c 

c 

cross-sectional area of the waveguide. 


is any 



Figure G.2. Waveguide geometry with fields radiated to the left (n>n + ) 
of S being of concern. 
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Although Equations (G.8) and (G.9) are given for a volume source 
distribution, it is very simple to modify them for a surface (or line) 
source distribution. For this purpose, the volume integral on the right 
hand side of Equations (G.8) or (G.9) should be replaced by a surface 
(or line) integral over the extent of the surface (or line) source 
distribution. 
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APPENDIX H 


EDGE CORRECTION FOR APERTURE INTEGRATION 


Z 



Figure H.l. Scattering from a perfectly-conducting half-plane. 

As shown in the figure, a perfectly-conducting half -plane is 
located in the region y=0, x<0. The geometry is infinite in the 
z-di recti on, therefore the problem can be reduced to a scalar problem in 

A 

terms of the z-component of the electric field (E z ) or the z-component 
of the magnetic field (H z ). The former case is called the soft case and 
the latter is called the hard case. The analysis for both cases are 
very similar, therefore only the soft case will be considered here. 
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Let S^, S^, denote the top and bottom surfaces of the half-plane, 
respectively and S denote the aperture surface making an angle 9 with 
the x-axis as shown in the figure. The half -plane is illuminated by a 
plane-wave 

E 1 . eJkpcosU-n (6 < < ,) (9 >|) (H.l) 

In (H.l) 0 and <t>' is restricted so that there exists reflected 
fields at the aperture surface Sa. The analysis can be simlarly carried 
out to the cases where this restriction is not present. The total field 
E z can be written as the combination of three contributions; namely, the 
incident, the reflected and the diffracted fields. The incident and 
reflected fields are also called the geometrical optics (GO) fields. 

The GO field is given by the following expression: 

E go . e jkpcos(*n , ). e jkpcos(»+V), e jk x x 2Js1nk ^ . 00 ^ ^ 

=0 ; on S b (H.2) 

since is in the shadow region of the half-plane. In (H.2), 

k] = kcos<j>' an( j k^ = ksin<f>' (H.3a;b) 

x y 

The transverse fields can be calculated using 

— a 3 a 3 

-jkZ Q H = (-y 97 + x 9p-) E z (H. 4 ) 

Therefore, the transverse geometrical optics fields can be obtained 
using Equations (H.2) and (H.3): 
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„ ,• jk n x 

-jkZ„ H g0 = jky e x 2cosk^y 


0 X 


(H.5) 




-jkZ Q H g0 = -jk^ e J x 2sink^y (H.6) 

In the physical optics (PO) approximation, one uses the radiation 
integral to get the scattered field using the GO fields on the 

illuminated side of the half-plane (S^). 

Therefore the PO approximation to scattered field is given by: 


E™ - / ds' [-jkZ Q G q (nxH go ) + (nxE go ) x V G Q ] 
S t 

* 9 * 3 . 

where VG q = (x gj + y g^J G q , and 


(H.7) 


G o (x,y | x' ,y' ) * "4 H q^ (Mp-p'I) 


(H.8) 


is the free-space Green's function in two-dimensions, and n is the unit 
normal to the surface. Using Equations (H.5) and (H.6) and the fact 
that E gu = 0 on S t , in (H.7), one obtains: 


E P0 = - / [- jkZ H • G ] dx’ 
Z Q J L 0 o x o J 


The Green's function G can also be written as: 

0 


6„ - - i Hk 


ly.y'l 

e J jk x' 

c e x 

x K y 


where 


k^ + k^ = k^ 
K x y K 


(H.9) 


(H.10) 

(H.ll) 
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Substituting (H.5) and (H.10) into (H.9) and changing the order of 


integration, one obtains: 

-jk x-jk v 
r P0 -1 f- e ■ r 

E z I dk x — r 


. » jkV + jk x 

k 1 / e x x dx' 

y o 


F P0 zl t 

E z = 2n J 


oo e 


-jk x x-jk y y 


k%k 

x x 


dk 


Using the transformation 


k x =kcosa 


k^=ksina , 


Equation (H.13) becomes 

E P0 _ zi r e -jkpcos(o-*) . lini ! — r d 

b z ~ 2n J e COSa+COS<(> act 


where the integration path is shown below: 



Rea 


Figure H.2. The integration path C . 


(H.12) 


(H.13) 


(H.14a;b) 


(H.15) 
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On the other hand, if one integrates the GO fields over S and S, to 

a b 

get the scattered field, one obtains: 

E* = J ds' [-JkZ ft G q (n x H go ) + (n x E go ) x VG q ] (H.16) 

a 


since GO field is zero on S. . 

D 

Using (H.2), (H.5), (H.6) and (H.10) in (H.16), one 
obtains: 

- j v j v ~ 

. J 00 e 

E* = - 4irsin9 J“ dk y Z [(- k v + k y ) sine+k cose] 

L -oo X a a y 

1 

k y + ^ k x +k x^ cot9+k x 




1 1 

- _ ; — 

k + (k 1 +k 1 cote+k 1 k + fk^+k 1 cote-k 1 
i_ y xx 7 x y ' x x' " y_ 


(H.17) 


Using the transformations in (H.14) and after some straightforward 
manipulations one finally obtains: 




-jkpcos(o-<|») — LLHi. — 
e cos a+cos <|> 


The Equations (H.15) and (H.18) are exactly the same, therefore the 
Ufimtsev edge correction [26] to (H.15) can be used without modification 
to correct for the aperture integration result of (H.18). 
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APPENDIX I 


APPLICATION OF THE RECIPROCITY THEOREM TO FIND THE p 
RELATION BETWEEN THE SCATTERING MATRICES [S 12 ] AND [S 2 i3 



Figure 1.1 The geometry of the problem. 


Figure 1.1 shows an open-ended waveguide cavity geometry with 

perfectly conducting walls. Let there be modal fields inside the 

* 

waveguide propagating in +n direction represented by: 

-+ + - "j S_n 

E + = C + (e . +e ) e n 
n n ' nt nn 


H + r + (Z - 
" = C n (h nt +h nn ) e 


( 1 . 1 ) 

( 1 . 2 ) 
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where e nt (or R ^) is the transverse electric (or magnetic) field, and 

I (or R ) is the longitudinal electric (or magnetic) field. C + is 
nn nn n 

the coefficient of the nth mode. 

Some of the energy carried by this modal field will be radiated into 
the free space through the aperture. The far field radiation can be 
calculated approximately by the aperture integration of the modal field 
and the resulting radiated fields are shown by E r and H r . 

The rest of the energy of the incident mode will reflect back to the 
guide in terms of an mth mode. 



(e -e ) 
' mt mn 




(1.3) 


H“ 

m 



(-R 

v mt mn 



Bn 

m 


(1.4) 


where r is the modal reflection coefficient from the open end. 

nm 

First consider the geometric plane defined by n=-L inside the 
semi -inf i nite waveguide and let S Q denote the area of this waveguide 
cross section at n=-L as shown in Figure 1.1. Then let S denote the 
surface area which tightly encapsulates the complete outer (exterior) 
surface of the semi -infinite waveguide, and also a portion of the inner 
walls of this waveguide up to the distance n=-L within the guide. The 
surface S does not include the plane at n=-L within the guide. Let E 
denote the sphere at infinity which surrounds the semi -infinite guide 
such that S on the exterior or outer wall of the semi -infinite guide is 
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connected to the surface Z at n+ -• as in Figure A.l. Next, consider 

the following two cases. In the first case, the semi -infinite 

rectangular waveguide is excited from within by the modal fields with 

amplitude C* which propagate in the +n direction. However, in the 

second case, the geometry is excited by external fields, E e and H e , 

which for convenience is assumed to be produced by an electric test 

source J at R exterior to and in the far zone of the semi -infinite 
e p 

guide. 

Let the equivalent sources J $ and M g be located at n=-L in the 
guide. 


J = n x 
s 


C + (R . +fi ) 
n ' nt nn 


je n L 

e + z r_ 
m 


mn 


C + (-R .+R ) 

n ' mt mn 




(1.5) 


= -n x 


~ . j B L . - j B L“ 

C + (e +e ) e n + z r C + (e . -e ) e m 
n nt nr mn n v mt mn 


( 1 . 6 ) 


where n is the unit vector pointing into the volume V, enclosed by 
surfaces S+S q +e. From reciprocity one obtains: 


ff E r x H e - E e x H r 

Z+S+S — 
o 


nds » /// 3 • E r dv 

V e 


(1.7) 


However, by the boundary conditions 


n x E' |=0 
on S 


( 1 . 8 ) 


and 


n x E c | = 0 
on S 


(1.9) 
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Also, (E r , H r ) and (E e , H e ) satisfy the radiation condition on E. 
Therefore (1.7) reduces to: 

- // H® • M ds + // E e • 5 ds = p . E r (R ) , (I. 10) 

S s S s e p 

o o 

since 

5 e - P e « ( |R-R p |) . (1.11) 

a point source in the far field. 

In (1. 11) Rp is the position vector from the reference point to 
point P. 

Inside the waveguide the fields E e and H e will have the following 
representation 



J3 n 


E e = E A" (e . - 
p P Pt 

■ 5 ) e p 
Po 

(1.12) 


•j 3 n 


H e = r A - i 

4 . R ^ O P 

n m 
\ *•*'*' / 

p "p v "pt 

"pn' 

Substituting (1.12) 

and (1.13) into ( I . 10) one obtains: 





-j 3 l . 

E A (-R +fi ) e P C + 
p P Pt po' n 


. j3_L _ . -J3 L“ 

e xn e + E r mn (e xn)e 
nt m mn mt 


ds 


~j 3 L + 

// E A (e.-e )e P C + 
S 0 P Pt Po' n 

o H 


J 3 L * - J* 3 L“ 

(nxR )e - E r (nxR )e 
nt m mn mt 

m 


ds 


- p e • E 


(1.14) 
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Using the orthogonality property of waveguide modes one obtains: 


V 'nt x h nt 
o 


n ds 


<- C n M> 2 = Pe • ^ V 


(1.15) 


P e • E < R p> 

■2 ct a" = : : — 

n n // e x h • n ds 
c n n 


(1.16) 


~ i A i ^ i 

In order to produce an incident field E = 0 E Q + 9 E^ a t the opening of 
the waveguide, one can let 


a -• a .• 4lt 

p e • - (e * ♦ e;> 7ST R p e p 


(1.17) 


where e and <f> are defined as the elevation and azimuth angles in a 
spherical coordinate system located at 0, in Figure 1.1. Thus, 


a ,• a , _ r - jkR 

(9 Eg ♦ ♦ e’) • E r (Rp) R p e 0 


A" = 


n -2 C + / / e . x R 

n rvh 


nt nt 


n ds 


4ir 

juu 


(1.18) 


In matrix notation, it is clear that (1.15) can be expressed as: 


[e£(P) eJ(P) e[(?)] 


ex 


ey 


ez 


* U^f) (-2[p nt xfi nt .Sds]) <[A-1) 

0 (1.19) 
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where (-2[ f/e nt .xR nt *nds]) 1s a dia 9 ona1 matrix, it is noted that 




E r contains C* since E r (P) is produced by the radiation of E*. 
Furthermore, it is obvious from (2.13) and (2.21) that 


E r (P) 
x ' 

E r (P) 

y 

eJ(p) 



( 1 . 20 ) 


and 



( 1 . 21 ) 


Incorporating (1.20) and (1.21) into (1.19) yields: 


i — — i 



p 

ex 


p 

ex 

([c;] t )([s 12 ]) t 

P 

ey 

= ([c;n T (-2[//5 nt xh- nt -Sds])as 2 P]) 

P 

ey 


P 

J 

0 

P 


ez 


ez 


( 1 . 22 ) 


The relationship between [Sj 2 !I and [S^] becomes evident from (1.22); 
namely 

[Si 2 ] T - (-2|/5 nt xfi nt .^mS 2 P] . (1.23) 

0 
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APPENDIX J 


EQUIVALENT MAGNETIC LINE DIPOLE 



Figure J.l. A wedge illuminated by a 8 Q -po1arized plane wave. 


In this appendix, the edge diffracted field from a wedge 

A. 

illuminated by a 3 0 -polarized plane wave (shown in Figure J.l) will be 
represented by the free-space radiation from magnetic line dipoles 
located at the edge position. The incident plane wave is in the 
following form: 


E 1 (s) = e 0 e 


a -i 

. ; .-ok 


(j.i) 
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where 


£ = -k(sine cos»j> x + sin 0 sin<j> y + cose z) (J.2) 

and i is the vector defining an observation point. The angles 0 1 and ^ 
are the elevation and the azimuth angles, respectively of the direction 
of the incident plane wave with respect to the wave with respect to the 
spherical coordinates. It is clear from Figure J.l that 0 1 =tt- 0 o . 

The electric field diffracted from the wedge and observed at the 
observation point P shown in Figure J.l is given by [26] 

A 

E d = -3 D s (i>,<f» 1 ;0 o ,0,a )\]j^) e -J ' ks E 1 (J.3) 

where D„ is the soft diffraction coefficient [ non-uni form soft 
s v 

diffraction coefficient is given in (3.3)), and p is the caustic 
distance for the diffracted rays. If p»s, then the diffracted field 
expression becomes 

2 . e' jks 

E d = -B D s (4>,4> 1 ;0 o ,0,ot) ~J=~ E 1 (J.4) 

Now, consider the free-space radiation from magnetic line dipoles 

A 

located along the z-axis as shown in Figure J.2. For an observation 
point with an elevation angle 0, the dipoles are given by 

M d = (x sin0 + z cos0 cos«|>) M d e "j kzcos0 fi(x)6(y) (J.5) 

where 0 is elevation angle and «|» 0 is a constant angle and S(x) is the 
Dirac delta function. 
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Figure J.2. Magnetic line dipoles along z-axis. 


The electric vector potential F is given by 

F = (ksine + z cose cos<j> 0 ) H^(k /x 2 +y 2 sine) e "^ kzcos0 


The radiated electric field is given by 

E = -vxF=$ {-kcosecos^sineH^+jkcosesinecos^H^^e - ^ 20050 

+ p yy { jkcosesi nesin<f> o H^ 2 ^ } e"^ 20050 
- z {ksin 2 esin<|» o H^} e" jkzcoS0 


where the arguments of the Hankel functions are as in (J.6). If one 
employs the large argument approximations to Hankel functions [44], 
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one obtains 


E = 




e-Jks 

/T 


(J.8) 


where 


s = /x 2 +y z +z^ 


(J.9) 


In (0.8) the angles 0 and <t> Q are replaced by 0 and <|>, respectively. 
Therefore, (J.8) applies for the observation point P of Figure J.l. By 
comparing (J.4) and (J.8), one deduces the strength of magnetic dipoles 
as follows: 


M 


d 


, /8ir D s ($,$"; 0 o ,0,o) ^ 

VW sTii» E (°) 


(J.10) 


289 



APPENDIX K 


DESCRIPTION OF THE EQUIVALENCE PRINCIPLE 
EMPLOYED IN THE CALCULATION OF SCATTERING MATRICES 


( E* ,H f ) 



Figure K.l. Radiation of a dipole dp in the presence of the 
semi -infinite waveguide. 

Let a dipole dp g radiate in the presence of the semi -infinite 

waveguide structure. The aperture of the waveguide is shown by and 

S denotes the surface of the waveguide wall. Let the total field be 
w 

represented by E t ,H t . By the equivalence principle one can use the 

equivalent sources J and M on S and S which radiate the total 

eq eq a w 

field inside the waveguide region and null field in the outside region, 
as shown in Figure K.2. 
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(E f , H f ) 
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Figure K.2. The equivalent problem of Figure K.l. 

The equivalent sources are given by 

5 eq = n x H 1 (K.l) 

H gq = E x n (K.2) 

A 

where n is the unit vector pointing into the waveguide region. By the 
boundary conditions, the equivalent magnetic sources on S is zero. If 
these sources are known, then the fields E t , H* coupled through the 
waveguide can be obtained. In order to calculate the fields one can 
react the sources and fields of Figure K.2 with the sources and the 
fields of Figure K.3, where a dipole dp & is radiating inside the 
infinite waveguide. That infinite waveguide is obtained by extending 
the semi -infinite wvaeguide of Figure K.l. 
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I 



Figure K.3. The dipole dp g radiating in the infinite waveguide. 

Clearly, the equivalent electric source 3 „ on S ., will not react 

eq w 

with the modal fields of the dipole dp g , since the tangential electric 
field in the modal expressions vanishes on S . Therefore, coupling 
through the waveguide modes will be determined by only the equivalent 
sources at the aperture S . Also, it is noted that since the equivalent 

a 

sources radiate null field external to the semi -infinite waveguide, the 
reaction of these fields with the dipole source of Figure K.3 will be 
zero. 
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APPENDIX L 


ILLUSTRATION OF RECIPROCITY IN THE APERTURE INTEGRATION 



Figure L.l. Modal radiation from the open end. 

In this appendix, the aperture integration technique will be used 
to obtain the modal radiation from and coupling through the open end of 
a semi -infinite waveguide in the front half -space. It will be shown 
that these two results satisfy the reciprocity relationship stated in 
Appendix I. 

The radiation problem is sketched in Figure L.l. The interior 
fields (E + ,H + ) propagate through the aperture at the open end and 
radiate into the exterior region. In order to simplify the analysis the 
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observation point is assumed to be in the far field of the aperture and 
separated by the distance r from the origin. The interior fields E + and 
H + are given in terms of a summation of waveguide modes, so 


z c 

n 



and 


(L.l) 


H + = Z C* H* 
n n n 


(L.2) 


where C* is the amplitude of the n^* 1 modal field. The aperture 
integration procedure assumes that the radiated field is due to the 


approximate equivalent electric and magnetic currents and located 
at the aperture. The currents and are given by 


J $ = n x H + 

(L.3) 

— - a 

M = E x n 

(L.4) 


where n is the unit normal vector at the surface of the aperture 
pointing into the exterior region. The approximate currents are then 
assumed to radiate in free-space. Therefore, the far-zone field 
radiated by the equivalent sources is written as 


;rad li 
E = 4 tt 


-jkr 


If 

S 

a 


jkr*r' 


[Z RxRxJ +RxM ] ds' 

OSS 


(L.5) 


where the unit vectors r and R are in r and R directions, respectively 
as shown in Figure L.l, and r' is the vector from the origin to the 
source point. The substitution of (L.l) through (L.4) into (L.5) yields 
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ik e 


i rad = 


•jkr 


E C 


" "I' 


jkr' *r 


[Z Q RxRx(nxH*)+Rx(E*xn)] ds' 


(L.6) 


The next step is to calculate the coupling of incident energy due 

to a small electric dipole dp & into the waveguide modes using the 

aperture integration. The dipole dp g is located at s 1 from the origin 

rn 

and the amplitudes of the modes excited by the dipole is shown by (C p ), 
The mode amplitudes are given by (3.97) 


u (i; • of - s; . ds' 

.GO _ h 

2 // E + x H + • ds' 

^ n n 


(L.7) 


where and are given by (3.92) and (3.93), 



Figure L.2. Coupling of incident dipole field. 
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which is the reciprocity relationship between the problems of the modal 
radiation and the coupling into the waveguide due to a dipole source. 
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